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^ ' Abstract 

D \ Ordinary-sterile neutrino oscillations can generate a significant lepton number asymmetry in 

'-^ \ the early Universe. We study this phenomenon in detail. We show that the dynamics of 

ordinary-sterile neutrino oscillations in the early Universe can be approximately described by 
^ ■ a single integro-differential equation which we derive from both the density matrix and Hamil- 

tonian formalisms. This equation reduces to a relatively simple ordinary first order differential 
equation if the system is sufficiently smooth (static limit). We study the conditions for which 
the static limit is an acceptable approximation. We also study the effect of the thermal dis- 
tribution of neutrino momenta on the generation of lepton number. We apply these results to 
show that it is possible to evade (by many orders of magnitude) the Big Bang Nucleosynthesis 
(BBN) bounds on the mixing parameters, and sin^2^0; describing ordinary-sterile neu- 
trino oscillations. We show that the large angle or maximal vacuum oscillation solution to the 
solar neutrino problem does not significantly modify BBN for most of the parameter space of 
interest, provided that the tau and/or mu neutrinos have masses greater than about 1 eV. We 
also show that the large angle or maximal ordinary-sterile neutrino oscillation solution to the 
atmospheric neutrino anomaly does not significantly modify BBN for a range of parameters. 



I. Introduction 



There are three main experimental indications that neutrinos have mass and oscillate. They 
are the solar neutrino problem||l[], the atmospheric neutrino anomaly^] and the Los Alamos 
LSND experiment . It is also possible that dark matter may be connected to neutrino 
masses[|]. The three experimental anomalies cannot all be explained with the three known 
neutrinos so it is possible that sterile neutrinos exist. 

A potential problem with any model which contains sterile neutrinos is that these extra 
states can contribute to the energy density of the early Universe and spoil the reasonably 
successful Big Bang Nucleosynthesis (BBN) predictions. For maximally mixed Ve and v'^ 
neutrinos and and v'^ (or and v'^) neutrinos (where the primes denote sterile species), 
the following rather stringent BBN bounds have been obtained p, ^ assuming that the 
lepton number asymmetry of the early Universe could be neglected: 

\Sml,\ ^ 10-« eV, \Sml^,\, \SmU ^ 10^' eV\ (1) 

Observe that if valid these bounds would rule out the large angle — u'^ oscillation solution 
to the atmospheric neutrino anomaly and would restrict much of the parameter space for the 
maximal oscillation solution of the solar neutrino problem|^, However, these bounds do 
not hold if there is an appreciable lepton asymmetry in the early Universe for temperatures 
between 1 — 30 MeV|TI|]. Remarkably, it turns out that ordinary-sterile neutrino oscillations 



can by themselves create an appreciable lepton number asymmetry [|T2 . 

The bound on the effective number of neutrinos N^^^ present during nucleosynthesis is 
the subject of some discussion recently. In Ref. |T^, it is argued that the current information 
suggests N^-l"-^ ~ 2.1 ± 0.3, while other authors dispute this conclusion. For example, in 
Ref . , Ref . [p!5| and Ref . |p!6| , the upper limits N^-^^ < 3.9, 4.5, 4.0 are respectively derived. 



Thus, it may be possible that N^-I'f = 4 is allowed. In this case note that many of the BBN 
bounds derived in Refs.0, ||, |^ |^, quoted in Eq.(0), need not apply. However, for the present 
paper we will assume that the bound on the effective number of neutrinos is less than 4. 
This is useful even if it turns out that A^^-^'^ > 4 is allowed. For example, the large angle (or 
maximal) ordinary sterile neutrino solutions to the atmospheric and solar neutrino problems 
may require N^-^-l" ~ 5 if they are to be solved simultaneously. Also note that in the special case 
of mirror neutrinos [p!7|] , the mirror interactions can potentially bring all three mirror neutrinos 



(as well as the mirror photon and electron positron pair) into equilibrium (equivalent to about 
6 additional neutrino species) if any one of the mirror neutrinos is brought into equilibrium 
above the neutrino kinetic decoupling temperature. 

The purpose of this paper is two-fold. First, we will study the phenomenon of lepton num- 
ber creation due to ordinary-sterile neutrino oscillations in more detail than in the previous 
studies |]T^, |T^. For example, we will study the effect of the thermal distribution of neutrino 
momenta. Using these results we will then study the issue of whether or not the generation 
of lepton number due to ordinary-sterile neutrino oscillations can reconcile the large angle 
ordinary-sterile neutrino oscillation solutions to the solar neutrino problem and atmospheric 
neutrino anomaly with BBN. 
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The outline of this paper is as follows. In section II, we discuss lepton number generation 
in the early Universe by ordinary-sterile neutrino oscillations and derive a simple equation 
describing the evolution of lepton number. We expand the analysis of Ref. [0 and discuss in 
detail the approximations behind this analysis. In section III, we will use the density matrix 
formalism to derive a more exact equation describing the rate of change of lepton number 
which is applicable even when the system is changing rapidly (e.g. at the resonance). In the 
appendix, we show how the same equation can be derived from the Hamiltonian formalism. 
Using this equation we derive the region of parameter space where the much simpler equation 
derived in section II is approximately valid. In section IV the thermal distribution of the 
neutrino momenta is considered. In section V we study the effect of non-negligible sterile 
neutrino number densities. We then apply these results to obtain the region of parameter space 
where large neutrino asymmetries are generated. We also determine the region of parameter 
space for which ordinary-sterile neutrino oscillations (with < and for |(5m^| ~ 10~^ eV'^) 
are consistent with BBN. Our work improves on previous studies]^, |^, |^, ||, because these 
studies were obtained without taking into account either the neutrino momentum distribution 
or the result that ordinary-sterile neutrino oscillations create lepton number. In section VI 
we first briefly review the large angle ordinary-sterile neutrino oscillation solution to the solar 
neutrino problem. We then show that the generation of lepton number due to ordinary-sterile 
neutrino oscillations can significantly relax the BBN bounds for this solution to the solar 
neutrino problem. We also show that the large angle or maximal z/^ — Us oscillation solution 
to the atmospheric neutrino anomaly is consistent with BBN for a range of parameters. In 
section VII we conclude. 

II. Lepton number creation from neutrino oscillations - static approximation 



Together with M. Thomson, we showed in Ref.|T^ that ordinary-sterile neutrino oscillations 



can create a large lepton asymmetry in the early Universe ||19||. A simple differential equation 
describing the evolution of lepton number was derived which seemed to work very well. We 
also checked our results with the more exact density matrix formalism ||2l|. Further numerical 
work, and analytical work based on the density matrix formalism, has subsequently been done 
in Ref. which confirms our results. 

For ordinary-sterile neutrino two state mixing, the weak-eigenstates [ua, i^s) will be linear 
combinations of two mass eigenstates {ua, ^b)'- 

= cos O^Va + siu = - siu O^Va + COS O^Vb. (2) 

Note we will always define 9q in such a way so that cos2^^o ^ (this can always be done). We 
also take the convention that 5m^^ = — m^. Hence with this convention 5m^^ is positive 
(negative) provided that rrif, > nia {nib < nia). 

In this section we will for simplicity neglect the effects of the thermal distribution of 
momentum, and assume that all of the neutrino momenta are the same and equal to the 
average momentum (i.e. p = (p) ~ 3.15T). In section IV we will consider the realistic case 
where the neutrino spread is given by the Fermi-Dirac distribution. Following Ref. fl^ , we 



2 



can derive a simple equation for the rate of change of lepton number due to colhsions and 
oscillations. Note that it is possible to identify two distinct contributions to the rate of change 
of lepton number. First, there are the oscillations between collisions which affect the lepton 
number of the Universe because neutrinos and anti-neutrinos oscillate with different matter 
oscillation lengths and matter mixing angles in the CP asymmetric background. Second, 
there are the collisions themselves which deplete and Ve at different rates. This is because 
the rates depend on the oscillation probability. The oscillation probability for ordinary-sterile 
neutrino oscillations is different to the oscillation probability for ordinary-sterile anti-neutrino 
oscillations (which is again due to the CP asymmetric background). Generally, the rate of 
change of lepton number is dominated by collisions in the region where the collision rate is 
larger than the expansion rate 1 12]. [A possible exception to this is in the resonance region 
where the matter mixing angle changes rapidly]. For the case of — fs oscillations (where 
a = e, /i, r), the rate of change of L^^ due to collisions is governed by the rate equation. 



dt 



+^r(z/. 



+ ^r(z7, ^ 



(3) 



where the n's are number densities and Ly^ = {ny^—nD^)/n^. Using T{ua J^s) = r('^s 
and r(i/Q — > Ug) = r(z/s — >• z/q,) (will justify this in a moment), Eq.(^) simplifies to 



dL,. 



dt 



n 



Vs 



(4) 



This equation can be re-written in the form 

dt 
where 



'^l - ^u) [-r(z/« ^ ^s) + r(z7, ^ z/,)] - (AC; - AC:) [r(^a ^ ^s) + r(z7. 



2n. ' 



2n. 



(5) 
(6) 



"7 "'"7 

Observe that ordinary-sterile neutrino oscillations do not change the total particle number, 
from which it follows that 

AC:+Ac-: = o. (7) 

Using Eqs.(|^|^, the rate of change of L^^ due to collisions is given by 



(8) 



where we have used rty^ +np^ ~ 3n^/4+0{Ll^). We will assume for the present that negligible 
sterile neutrinos are produced, i.e. ny^,np^ -C ny^,np^, and hence AC" <S 1- 



In order to work out the reaction rates, we can invoke a simple physical picture] 22, 23 



24| . The oscillations of the neutrino between collisions produce a superposition of states. 
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The collisions are assumed to collapse the wavefunction into either a pure weak eigenstate 
neutrino or a pure sterile eigenstate neutrino. In other words, we assume that the collisions 
are measurements (in the quantum mechanical sense) of whether the state is a sterile or weak 
eigenstate. The rate of the measurements is expected to be the collision frequency T^^. 

Actually it happens that the above picture is not completely correct. It turns out that 
it does lead to an accurate description only if the rate of measurement is taken to be half of 
the collision frequency that a pure Ua state would experience [p^ ]. This applies to both sterile 



neutrinos and ordinary neutrinos. Thus using this result the reaction rate T{i/a i^s) is given 
by half the interaction rate of the neutrino due to collisions with the background particles 
multiplied by the probability that the neutrino collapses to the sterile eigenstate that is 

rK-^.) = ^(P..-...). (9) 

The thermally averaged collision frequencies T^^ are 

- VaGlT^ (10) 

where ye ~ 4.0,?/^,^ ~ 2.90, Gp is the Fermi constant {Op ~ 1.17 x 10^^^ MeV'^) and T 
is the temperature of the Universe [equations analogous to Eqs.(^|rU[) hold for antineutrinos]. 
The quantity Py^^y^ is the probability that the neutrino Va collapses to the sterile state Vg 
after a measurement is made. The brackets (...) denote the average over all measurement 
times. Note that Pu^^us = Pua^v^^ so it follows that T{va ^s) = ^{^s ^a) (given that 
the rate of measurement is the same for ordinary and sterile neutrinos p^) and similarly for 
the anti- neutrino rates. In the adiabatic limit, 

(P,„^.J^sin2 2^„(sin2-^), (11) 

where r is the distance (or time) between collisions. The quantities 6m and Lm are the matter 
mixing angle and matter oscillation length respectively. They are related to the vacuum 



parameters Oq and Lq by l23, 26 



9 ^ sin^ 29n 

8111^26^ = ^ ^, 12 

l-2zcos29o + z^' ^ ' 

and ^ 

^"^ ^ VI -2z cos 2^0+^' ^^^^ 

where l/i^o = Aq = 5m? I2p. In this equation, z = 2pVa/Sm'^ where Va is the effective 
potential due to the interactions of the neutrinos with matter and p is the neutrino momentum. 
The effective potential is given by 

Va = i-aP + If)Al (14) 
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where the dimensionless variables and are given by 

Al AgM^ {py 

where Mw is the W-boson mass and Ae ~ 55.0, A^^r — 15.3 (note that the "p" superscript 
serves as a reminder that these quantities are neutrino momentum dependent). The function 
is given by 

= L,^ + L,^ + L,^ + L,^ + 77, (16) 

where ?7 is a small asymmetry term which arises from the asymmetries of baryons and electrons. 
It is given by|25] 

= (^ + 2 sin^ d^)Le + - 2 sin^ d^)Lp - ]^Ln ^ ^Ln, (17) 

where sin^ 6'^, is the weak mixing angle and we have used = Lp ~ L^r. Thus r] is ex- 
pected to be of order 10~^°. Note that the matter mixing angle 9m and oscillation length Lm 
for antineutrino oscillations are obtained from Eqs.( p!2| - p!5| ) by performing the transformation 



We denote the thermal average of the variables a^, If hj a = (a^), b = (If). From Eq. ([l5|) , 
they are given approximately by 



5m? \ 5m? j \ MeV I ' 

n2 / N 6 



^ ~G.ZV2TGpn A^T^ _ _ / T ( eV±\ r _ QSCiUationS 

5m? Ml, — Vl3 Mev) \5m^ ) ^ ^'^^ '^^ '^^ OSCUiailOnS, 

-6.3V2TG^n A,,.r^ _ _ / T ^ ( sYl] for z/ - u oscillatious 

<5m2M^ — \i6 MeV ) U^-V' ^'^ z/^ oscmaxions. 



where we have used = 2({3)T^ /tt"^ ~ T^/4.1 [C(3) — 1.202 is the Riemann zeta function of 
3]. The matter mixing angles 6^, dm expressed in terms of the parameters a, b are given by 

sin^ 26^ = T^—jT sin^ 29^ = . ^ ^ (u ^ ^' (^9) 

[s^ + [b — a — cY\ [s^ + (o + a — cY\ 

where s = sin26'o, c = cos 2^0 • A resonance occurs for neutrinos when 6^ = 7r/4 and 
for antineutrinos when 6^ = 7r/4, which from Eq.(|19D implies that b — a = cos2^o and 
b + a = cos 26q respectively. In our analysis we will often need to consider the two distinct 
cases of very small mixing and very large mixing. For small mixing, cos26'o — 1 and the 
resonance conditions become 6 — a ~ 1 and 6 + a ~ 1. For large mixing, cos 2^0 — and the 
resonance conditions become a ~ 6 and —a ~ b. 

Using the above analysis, we can derive a simple equation for the rate of change of L,^^, 



I'J-'Ua _3_ 

dt ~ 16 



r„ 



- sin2 2ft^(sin2 (^)) + sin^ 20„(sin2 [^)) 
sin^ 2em{sm^ {^) ) + sin^ 2^^(sin2 (^) )] . (20) 
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The function, (sin^ (2!")) given by 

where Uq = 2tq = is twice the mean time between colhsions (of a pure weak eigenstate) 

and t is the age of the Universe (note that t ~ 00 is a good approximation because u^o ^ t)- 
Evaluating Eg. (pi]) we find 




sin^ h^ ) = - , L (22) 



where we have assumed that uq and Lm are approximately independent of t (static approxi- 
mation). Thus, using Eqs.(^, |^), we can rewrite Eq.(pOD in the form 



dt 8[x + {c-b + a)^][x + {c-b-ay] 
where A is a small correction term 



+ A, (23) 



. ^ 1 L.,5^r.Jx+(c-&)^ + a^] 

2[x + {c-b + a)^][x+{c-b-ayy ^ ' 



^2 I 7^2 / -f^ I „2 I o in-iy I 1 



and a; is given by 

where we have assumed p = (p) ~ 3.15T in deriving the last part of the above equation. 
Note that the correction term Eq.(P^) is smaller than the main term [Eq.(^)] provided that 
\LyJ\ <^ \a\. In the region where the correction term is larger than the main term, its effect 
is to reduce \LyJ\ such that Ly^ 0. From Eq. ([T8|) , the condition \LyJ\ > \a\ only occurs for 
quite low temperatures, 

T < 1 /|(5m2i^ ^ 



MeV 3 V eV^ 



(26) 



From the above equation, we see that in the main region of interest (T ~ 3 MeV), the 
correction term is much smaller than the main term provided that |(5m^| ~ 10^ eV'^. Note 
that for very large \6m^\ ~ 10^ eV'^, the correction term may be important. 

Observe that Eq.(PBD differs slightly from the equation derived in Ref.|T2|. The difference 
is that in Ref.[|^, we assumed that cUq/L^ ^ 1 (so that (sin^ r/2Lm) ~ 1/2) which is always 
true except possibly at the very center of the resonance |T2]. Also note that in Ref. [T^ 
we neglected a factor of 2 which arises because we negligently assumed that the rate of 
measurement was equal to the rate of collision. 

We now pause to review and comment on the assumptions made in deriving Eq. ([23|) . 
There are five main simplifying assumptions: 
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(1) We have neglected the thermal spread of the neutrino momenta, and have replaced all 
momenta by their thermal average (p) ~ 3.15T. 

(2) We have assumed that nu^,np^ <^ ny^,np^. If the number densities n^^^Uy^ are non- 
negligible, then we must multiply the first term on the right-hand side of Eq.(PB|) by the 
factor [n^^ - n^J/n^,^. 

(3) We have assumed that the transformation from the vacuum parameters to matter param- 
eters i.e., sin^o ~^ sin^^ and Lq — > diagonalizes the Hamiltonian. This is only strictly 
true in the adiabatic limit {\d9mldt\ -C |Am|). In the general case p6| . 



(2J 



with 

dOm 1 sin 29q d(h — a) 

"dT ~ 2(6-a-cos2^o)^ + sin2 26^ Jt ' 

where v}^ are the instantaneous matter eigenstates and = ^/L^- Expanding out 7 = 
\{d6m/ dt) / we find (neglecting da/dt), 

7 < 2(4) X 10-8 (^i^)'^' (pSi)'^' ' away from resonance, 
7 ~ 2(4) ( gj^° 200 ) (pm^) ^ ' initial resonance where b = cos 26^0, a ~ 0, 

7 ^ 6(9) X 10-^r=^ {^^) (|f^) at the resonance where \b - a\ = cos2^o, (29) 

for z/g — {v^^T ~ ^s) oscillations. However at the initial resonance where h = cos2^0) ct — 0, 
Liy^ is created rapidly. The contribution to 7 from a rapidly changing L,y^ at this resonance is 



2 



^ sin2 2eo \\Sm^\J d{T/MeV)' ^ ^ 

for z/g — i^^i,T ~ ^s) oscillations (and we have assumed that cos 26*0 ~ 1). Thus away from 
the resonance the adiabatic approximation is valid for the parameter space of interest, (i.e. 
for l^m^l ~ 10"^ eV^). However at the resonance the adiabatic approximation may not be 
valid. 

(4) Equation (p3D neglects flavour conversion of neutrinos passing through the resonance (the 
MSW effect). Observe that there is not expected to be significant flavour conversion at the 
initial resonance (where b ~ cos2^o) due to the MSW effect (even if the system is adiabatic 
at this resonance) because the frequency of the collisions is such that (sin^r/2Lm) <^ 1 at 
the center of the initial resonance, for most of the parameter space of interest. Indeed, at the 
center of the resonance, 

t^O _ sin 200 (5m^ ~ A v 1 fl^ cln 9fl„ f MsX^ ^ ^"i^ 

2L„ - ycGj,T5 2p - ^ J-U sm Zt7o y J , 

~ 90 tan 2^0, if & = cos 2^o- (31) 
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Thus, for sin^ 2^0 ~ 10~^, (sin^r/2Lm) -C 1. Note however that for temperatures below the 
initial resonance, the MSW effect may be important if there are neutrinos passing through 
the resonance. 

(5) We have assumed that the rate of change of lepton number is dominated by collisions. 
There is also a contribution from oscillations between collisions. Oscillations between collisions 
affect lepton number because the oscillations produce a superposition of states, where the 
averaged expectation value of the state being a weak-eigenstate is 1 — sin^ 26^. for neutrinos. 
This probability is generally unequal to the analogous quantity for anti-neutrinos, which is 
1 — sin^2^m- It is possible to show||12|| that for temperatures greater than a few MeV, the 
change in lepton number due to the oscillations between collisions is generally smaller than the 
change due to collisions except possibly at the resonance where sin^ 29m is changing rapidly. 

The effect of the thermal spread of the neutrino momenta should be to make the creation 
and destruction of lepton number much smoother. At any given time, only a small fraction of 
the neutrinos will be at resonance (because the resonance width is much less than the spread 
of neutrino momenta). Thus, the regions away from resonance may also be important. We 
will study the effect of the thermal distribution of momenta in section IV. 

The second assumption [(2) above] will be approximately valid for much of the parameter 
space of interest. This is because we are essentially interested in the region of parameter 
space where the sterile neutrinos do not come into equilibrium with the ordinary neutrinos. 
We will study the effect of the sterile neutrino number density being non-zero in section V. 
The assumptions (3) and (5), may not be valid in the resonance region. Note that we will 
denote the assumptions (3) and (5) collectively as the static approximation because in limit 
where the system is sufficiently smooth they will be valid. 

Clearly a more exact treatment of the resonance is desirable, since assumptions (3) and 
(5) may not be valid there. In section III we will develop a more exact treatment of the 
resonance region by examining the appropriate equations from the density matrix. As we will 
show in section III, this treatment leads to the following equation for the rate of change of 
lepton number 



dt ~ S 



sm 



\-dt" 



dr. 



where 



6rn? 



sin2^n, 



[cos2^o -&), A" 



-a. 



(32) 



(33) 



2p ' 2p ' ' 2p 

This equation is valid given the assumptions (1), (2) and (4) but does not require assumptions 
(3) and (5) [above]. This equation is an integro-differential equation and although compact 
cannot be solved analytically except in various limits. Note that the static limit corresponds 
to taking A^ as constant (i.e. independent of t'). In this limit Eq.(|32D reduces approximately 
to Eq.(^) as expected. In the appendix we show that Eq.(P^ can also be obtained using the 
Hamiltonian formalism provided that the rate of measurement is taken to be half the collision 
frequency. 
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Qualitatively, it turns out that the simplified equation, Eg. , gives a reasonable de- 
scription of the creation of lepton number as the Universe evolves. Assuming that Eq.p^) 
is valid, we now analyse the behaviour of L,^^ as driven by — Vs oscillations in isolation. 
Suppose that all initial asymmetries other than Ly^ can be neglected so that L^"-* ~ 2Lj,^. 
Notice first of all that for > if follows from Eq.(^) that h is negative and a has the 
opposite sign to Lj^^. Thus from Eq.(^) it is easy to see that the point L,^^ = is always a 
stable fixed point. That is, when L^^ > the rate of change dL^^/dt is negative, while when 
Ly^ < the rate of change dL^^/dt is positive, so Ly^ always tends to zero. [In the realistic 
case, where the baryon and electron asymmetries are not neglected, L'-"^ is given by Eq. ([l6|) . 
In this case is an approximate fixed point. Note that even if all of the lepton numbers 

where initially zero, lepton number would be generated such that L*^") ~ 0, i.e. 2Ly^ ~ —77 
[see Eq.(|TBp]. Note that 2Ly^ is only approximately —77 because of the A term proportional 
to in Eq.(ll)]. 

Now consider neutrino oscillations with 5m? < 0. In this case h is positive and a has 
the same sign as Ly^. From Eq. (^3]) , L^^ ~ is a stable fixed point only when h > cos2^o- 
When b < cos 26^0, the point Ly^ ~ is unstable. [That is, if L^^ > 0, then dLy^/dt > 0, 
while if Ly^ < then dLy^/dt < 0]. Since b ~ T^, at some point during the evolution of the 
Universe b becomes less than cos26'o and Ly^ = becomes unstable. If \5m^\ ^ 10"^ eV^, 
then this point (where b = cos2^o) occurs for temperatures greater than about three MeV 
(assuming cos2^o — !)• In this region the rate of change of lepton number is dominated 



by collisions and Eq.(p3D is approximately valid. When the critical point where b = cos2^o 
is reached, the lepton asymmetries are small and hence \a\ -C cos2^o — 1- Equation (^) 
then implies that dLy^/dt is approximately proportional to Ly^, which leads to a brief but 
extremely rapid period of exponential growth of Ly^ ||T^. Furthermore note that the constant 
of proportionality is enhanced by resonances for both neutrinos and antineutrinos at this 
critical point (a ^ 0, 6 = cos 26*0). The exponent governing the exponential increase in Ly^ 
is thus a large number (unless sin^ 29o is very small). Note that the critical point b = cos26'o 
occurs when 



T.. 13(16) fM!^^^^) MeV, (34) 



for the Ue — Ug {i'^i,t — ^s) oscillations we have been focusing on. 

As the system passes through this critical temperature, lepton number is rapidly created 
until a ~ cos 2^0 ~ The resonance at a = cos2^o ~ b cicts like a barrier which keeps 
a > cos 26q — 6 as the temperature falls below Tc. Since the parameter a is proportional to 
Ly^T'^, it follows that the lepton number continues to grow approximately like after the 
resonance as the temperature falls. 

As the temperature drops, eventually the oscillations cannot keep up with the expansion of 
the Universe. For temperatures well below the resonance, a ~ cos 26q (assuming that Ly^ > 
for definiteness) . In this region, the rate of change of a due to the oscillations is balanced by 
the rate of change of a due to the expansion of the Universe. That is, 

da _ _d^dLy^ da ^ 

dt ~ dLy^ dt ^ dt- ^^^^ 
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Eventually, the rate of change of a due to the expansion of the Universe becomes larger in 
magnitude than the maximum rate of change of a due to oscillations. At this point, a falls 
below the resonance point (i.e. a < cos2^o ~ b) and the value of L^^ will be approximately 
frozen. The point in time when this occurs is thus governed by the equation 

da dLy„, , da 

(36) 



dL,^ dt """" dt' 

The maximum rate of change of Li,^ occurs at the resonance where a = cos 26*0 — b. Using 
Eq.(p3l), we can easily evaluate dL^^/dt at this point. Assuming that cos26'o — 1, we find at 
the resonance, 

= — r^^a, 37 
dt 32 ^ ' ^ ' 

where we have assumed that x ~ sin^ 26q, which should be valid since we are in the region of 
low temperatures T ~ 3 MeV [recall that x is defined in Eq.(p5D]. Also note that 



da dadT iad.bT^ 



dt dT dt T Mf 



(38) 



where we have used the result that the parameter a is proportional to T^, and dT/dt ~ 
-5.5T^/Mp (which is approximately vahd for 1 MeV ~ T ~ 100 MeV, and Mp ~ 1.2 x 
10^^ MeV is the Planck mass). Thus, using Eqs.(^,^), the condition Eq.(^) can be solved 
for T. Doing this exercise, and denoting this value of T by Tf, we find 



Tf 



505m^ 


7 


5m? 









MeV, (39) 



where we have used Eq. ([To|) , Eq.(^). Thus, we expect Ly^ to evolve like until quite low 
temperatures of order 1 MeV. Note however that when the momentum distribution is taken 
into account, the situation is somewhat different. This is because only a small fraction of 
neutrinos (typically of order 1 percent or less) will be at the resonance, so that the magnitude 
of the maximum value of dL^^/dt will be reduced by a few orders of magnitude. Because of 
the 1/7 power in Eq.(^), the temperature where L^,^ is approximately frozen, Tf, increases 
by only a relatively small factor of 2 or 3. Finally recall that for temperatures below the 
initial resonance, the MSW effect can also contribute significantly. This is because for low 
temperatures near Tf, there will be a significant number of neutrinos which will be passing 
through the resonance. For low temperatures, the adiabatic condition is expected to hold [for 
most of the parameter space of interest, see Eq.(pPD]. Also, recall that the oscillations will not 
be damped by collisions for low temperatures [see Eq.(^)] and thus ordinary neutrinos can 
be converted into sterile neutrino states simply by passing through the resonance[^]. This 
effect will help keep a ~ 1 for even lower temperatures. 

Clearly these factors (the momentum distribution and the MSW flavour conversion of the 
neutrinos passing through the resonance) will be important if one wants to know the final 
magnitude of L^y^. For example, the final magnitude of L^,^ is very important if one wants 
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to calculate the region of parameter space where the L^,^ is large enough to affect big bang 
nucleosynthesis through nuclear reaction rates. However, for the application in this paper, 
the precise value of L^^ at low temperatures is not required, so we will leave a study of this 
issue to the future. 

In order to illustrate the evolution of L,y^ we take some examples. It is illuminating 
to compare the evolution expected from the simple Eg. p3|) [based on the assumptions (1)- 
(5) discussed above], with the evolution governed by the more complicated density matrix 
equations. [Eqns.(P6|), see next section for some discussion of the density matrix formalism]. 
The evolution of L,y^ as governed by the density matrix equations hold more generally than 
Eq. (p3D . This is because they do not require the assumptions (2), (3), (4) or (5) [discussed 
above] to hold. They do still incorporate assumption (1), that is the thermal distribution of 
the neutrino momentum is neglected. 

In Figure 1,2 we plot the evolution of L,^^ for some typical parameters. We consider for 
example t — Ug oscillations. In Figure 1 we take (5m^ = —1 eV"^, and sin^ 26q = 10~^, 10~^. 
Figure 2 is the same as figure 1 except that = —1000 eV^ and sin^ 26*0 = 10~^, 10~^. The 
solid lines are the result of numerically integrating the density matrix equations, while the 
dashed lines are the results of numerically integrating Eq.(p3D. We stress that in both the 
density matrix equations and in Eq.(^), the momentum distribution of the neutrino has been 
neglected. The effect of the momentum distribution will be considered in detail in sections 
IV, V. 

In the examples in Figure 1,2 the initial lepton asymmetry was taken as zero. The gen- 
eration of lepton number is essentially independent of the initial lepton number asymmetry 
provided that it is less than about 10~^[En, ITT|. This is because for temperatures greater than 



the resonance temperature, the oscillations destroy or create lepton number until L^"^ ^ 
independently of the initial value of L^^ [which we denote as Linu], provided that \Linit\ is 
less than about 10~^. For \Linit\ ~ 10~^, the oscillations at temperatures above the resonance 
temperature are not strong enough to destroy the initial asymmetry. Consequently, L^^ re- 
mains large, and it will become larger due to the oscillations which create lepton number at 
temperatures below the resonance temperature. 

As the Figures show, the behaviour expected from Eq.(p3D occurs. The main difference 
arises at the resonance where the magnitude of the lepton number is somewhat larger than 
expected from Eq. (p3|) . This occurs because the assumptions (3) and (5) [discussed above], 
which lead to Eq.(^) are not valid at this resonance. Actually, in Figure 1,2 we have plotted 
|. Integration of the density matrix equation reveals that in example 1 (but not in example 
2, although L^^ does change sign), the generated lepton number oscillates at the resonance 



and changes sign a few times (see Ref.|]T8[, for a figure illustrating this). Note that this effect 
can be understood from Eq . (p2D . To see this, observe that when Lj,^ is initially created at 
the resonance, the parameter A~ grows very rapidly because it is proportional to L,y^. The 
creation of L,^^ may be so rapid that //_^ X^dt' is approximately independent of r when the 



initial rapid growth of L^^ occurs. If this happens then at this instant Eq.(p2D can be simplified 
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to the approximate form 



dt ^ 8 



sm 



X-dt' 



t—WQ 



\ + dt" 



dr. 



(40) 



The oscillations occur because of the factor sin X'dt' which oscillates between ±1. Note 
however, that this oscillation of lepton number would not be expected to occur in the realistic 
case where the thermal spread of neutrino momenta is considered. 

Note that it may be possible to predict the sign of the asymmetry in principle. Assuming 
that the resonance is smooth enough so that Eq.(^) is valid, the equation governing the 
evolution of Ly^ has the approximate form 



dL^ 
dt 



(41) 



where V= rj + L^^ + L^^ + — L^^ (we have defined V such that it is independent of L^^). 
Note that A and B [which can be obtained from Eq. (|23|)] are complicated functions of time. 
Observe however that -B > and A is initially less than zero, and at the resonance A changes 
sign and becomes positive after that. In the region where 2A < B, the lepton number evolves 
such that 

{2A - B)L,^ + AV^O. (42) 

Thus will evolve such that it has a sign opposite to V just before the resonance. When 
2A > B, Ly^ will become unstable and grow rapidly. Note that at the point A = B/2, 



dL, 



dt 



A V . 



(43) 



Hence, at the point where the initial rapid creation of L^,^ occurs, the rate of change of L^^ 
will be proportional to V- Thus, we might expect that the sign of L^^ will be the same as 
the sign of the asymmetry V after the resonance. This means that L^^ should change sign 
at the resonance. Note however that because V depends on the initial values of the lepton 
asymmetries which are unknown at the moment, it seems that the sign of L^^ cannot yet be 
predicted. However the above calculation shows that the sign of L^^ should not depend on 



statistical fluctuations, as we initially thought likely 1 12 



Finally we would like to comment on the region of parameter space where significant 
generation of lepton number occurs. Firstly, we require that < and that |5m^| ~ 

10"^ eV^, so that ~ 3 MeV. For \6m^\ ~ 10"^ eV^, lepton number can still be generated 
but it is dominated by the oscillations between collisions and is oscillatory ||18|, |28| . Note that 
in the realistic case where the spread of momenta is taken into account, oscillations of lepton 
number would be smoothed out and may not occur. A numerical study in Ref.|jT8| shows 



that sin^2^o ~ 10 ^^(el^^/|5m^|)^/^ is also necessary (see also Ref. |]T2| for an approximate 



analytical derivation). Finally, we must require that sin^2^o be small enough so that the 
sterile neutrinos do not come into equilibrium. [For example, if there are equal numbers of z/^ 
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and z/g then the rate n^^T{i'^ ~^ K) = '^v'J^iy'e ~^ ^m) -^m cannot be generated]. We will 
re-examine the region of parameter space where significant generation of lepton number occurs 
in section V (where the effects of the Fermi-Dirac momentum distribution of the neutrino will 
be taken into account). 



Note that in Ref . [p^ , it is argued that lepton number generation only occurs provided 
that |(5m^| ~ 100 eV"^. We have not been able to verify this result, either analytically or 
numerically. In fact, we have been able to obtain any significant upper bound on |5m 
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III. Lepton number generation due to neutrino oscillations - A more exact treat- 
ment 

In this section we derive a more general equation describing lepton number generation in 
the early Universe which can be applied when the system is changing rapidly, as occurs, for 
instance, at the resonance. The only assumptions that we will make are the assumptions (1), 
(2) and (4) (discussed in the previous section). That is we will neglect the spread of neutrino 
momenta and set (p) ^ 3.15T, and we will also assume that there are negligible numbers 
of sterile neutrinos generated. In the appendix an alternative derivation (with the same end 
result) based on the Hamiltonian formalism is presented. Although not yet realistic because 
of assumptions (1) and (2), this derivation turns out to be particularly useful because it allows 
us to work out the region of parameter space where the simple Eq. ( |23D is approximately valid. 
As we will show, it turns out that Eg. (^31) has a wider applicability than might be expected 
from the adiabatic condition Eqs.(^Dp. 

The system of an active neutrino oscillating with a sterile neutrino can be described by 
a density matrix. See, for example, Ref. PT| for details. Below we very briefly outline this 
formalism and show how it leads to an integro-differential equation which reduces to Eq.(|2^) 
in the static limit. 

The density matrices for the neutrino system are given by 

Pu = ^ , Pu = ^ , (44) 

where Pq and Pq are the relative number densities of the mixed neutrino and anti-neutrino 
species, and P and P are the polarization vectors that describe the internal quantum state of 
the mixed neutrinos in terms of an expansion in the Pauli matrices a. The number densities 
of z/q, and Ug are given by 

with analogous equations for the anti-neutrinos. The evolution of Pq, P are governed by the 



equations [21 



|P = V X P + (1 - P^)(|/nPo)z - {D^ + D' + fjnPo){P,5i + P,y) 

ftPo = Ei=e,z.^;/3^a(r(z/aZ>a ^ n))iKnana " n^.UpJ, (46) 
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where A;^ = 1 and Ag = 1/4, and (...) indicates the average over the momentum distributions. 
The quantity V is given by 

V = /?x + Az, (47) 

where /?, A are defined by 

13 = sin 2^0, A = (cos2^o-&±a), (48) 

2p 2p 

where the +(— ) sign corresponds to neutrino (anti- neutrino) oscillations. The quantities 
and are quantum damping parameters resulting from elastic and inelastic processes 
respectively. According to ref. + = The function (^{(pilj 0V)) is the 

collision rate for the process (pip — > p)'ip' averaged over the distribution of collision parameters 
at the temperature T assuming that all species are in equilibrium. 
Expanding out Eq. (|46|) , we have: 

^ = /3P, + (l-P,)(|/o^Po), 

^^ = AP,.-/3P,-P,/^o, 

^ = -\Py - Pjujo, (49) 

where uq = 1/{D^ + + ^logPo) ~ 1/-D (where D = + D^). If we make the approxi- 
mation of setting all of the number densities to their equilibrium values, and also assume that 
the number of sterile species is small, then P;^ ~ 1 and Eq.(H3) simplifies to 



dPz pp 

XP.,-/3-Py/u;o, 

XPy - Pjujo. (50) 



dt — ''--^ X I--" ^ y 

dt ~ 



Strictly speaking, the approximation of setting P^ = 1 = constant can only be valid when 
(3Pz is small enough, so that MSW fiavour conversion cannot occur, i.e. when 



\P\ ~ |A| or — . (51) 

It is useful to introduce the complex variable P (t) defined by P= P^ + iPy. It is easy to see 
that the resulting equation describing the evolution of P (t) is given by 

dp ^ p 

t-^ = -\P-i— + p. (52) 
dt ujq 

The solution to this equation with initial condition P (0) = is: 

P (t) = -% /*/3(t')e^*'"*^/"°e'/t'^''*"dt', (53) 
h 
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where cuq has been assumed to be independent of time which is approximately vahd for tem- 
peratures above a few MeV where the expansion rate is less than the collision rate. One can 
easily verify that Eq.(|53|) is indeed the solution of Eq.(p^ by direct substitution. Thus, taking 
the imaginary part of both sides of Eq.(p3), we find that 



- / /3e 



{t'-t)/u;o 



COS 



Xdt" 



dt'. 



From Eq. p5|) , it follows that: 



dt 



IQdt^^ 



(54) 



(55) 



where denotes the z-component of the polarization vector for anti-neutrinos. Thus using 
Eq.(pOD the above equation becomes 



dL, 



dt 



— {Py-Py). 

16^ ^ ^' 



(56) 



Note that Py is given by Eq. ( ^Ij) and Py is defined similarly to Py except that we must replace 
a —>■ —a. Thus, we obtain 



dL,, 



dt 



-3^ 
16 



COS 



Xdt" 



COS 



Xdt" 



dt', 



(57) 



where A = 5m^(c — b + a)/2p, X = 5rn?{c — b — a)/2p. Note that we have taken P outside 
the integral, which is valid for T ~ 2 MeV, because /3 is approximately constant over the 
interaction time scale t — Changing variables from t' to the variable r where t = t — t' , 

this equation reduces to 



dL,. 



-3/32 



dt 



16 Jo 



COS 



Xdt' 



t~T 



COS 



Xdt' 



t-T 



dr. 



or equivalently. 



dL^ 
dt 
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(cos / Xdt') — (cos / Xdt') 



t-T 



t-T 



Note that the above equation can be re-written using a trigonometric identity, so that 



dL,, 



dt 



e"^/"'" sin 



X+dt' 



t-T 



sm 



X-'dt" 



t-T 



dr. 



(5J 



(59) 



(60) 



where A± = (A ± A)/2. 

The phenomenon of neutrino oscillations can also be described by the Hamiltonian for- 
malism. We show in the appendix that this formalism also leads to Eq.(|60D under the same 
assumptions. The density matrix formalism is particularly useful if one wants to keep track of 
the various number densities. In this more general case, it is very difficult to solve the system 
analytically and so far this more general case has only been studied numerically. 
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In the static limit where A, A are approximately constant, it is straightforward to show 
that Eq.(|58D reduces to Eq.(p3[) with x given by 



1, 



rather than by Eq.(|25|) [note that Eq.(|6TD reduces to Eq.(|25|) for most of the parameter space 
of interest except for quite low temperatures]. This difference is due to the fact that in 
deriving Eq.(^S|) we have made the assumption Eq. (|^) . Because Eq.(^) is much simpler 
than Eq. (|60|) , it is particularly useful to determine the region of parameter space where the 
static hmit [Eq.(p3D] is an acceptable approximation. We now study this issue. 

Expand Xf (note that we are using the notation that A^ denotes A evaluated at the point 
x) in a Taylor series around the point t' = t, that is 

\,, = \, + [t'-t](^\] +... (62) 



Using this Taylor series, the integrals //„^ Xdt' can be expanded as follows (with a similar 
expansion for //„^ Xdt'), 



(63) 



The static approximation will be valid provided that 



(cos /* Xdt') - (cos /* Xdt') ^ (cos r At) - (cos r At). (64) 

Jt-T Jt~T 



Using the expansion Eq. (|63D , observe that 

T f dX' 



cos / Xdt') = (cost 

Jt-T 



"'-2 5? . + 



>. (65) 



The above equation can be used to determine the region of validity of the static approximation 
Eq. (|6^ . The region of validity of Eq.(|6^) depends on the values of the parameters A, A. There 
are essentially four regions to consider. 



(a) ciJolAtl, cJojAtl ~ 1. In this region. Eg . (|64D is approximately vahd provided that 



(b) tuolAtl ~ 0, uJo\Xt\ ~ 1. In this region, Eq.(|6^) is approximately valid provided that 
Eq.(|66l) holds and 

Xdt') ~ 0. (67) 



(cos 



t-T 



16 



From Eg. this equation implies that 




(6J 



(c) ciJolAtl ~ 0, ido\Xt\ ~ 1. In this region, Eq.(|6^) is approximately valid provided that Eq.(l6^) 
holds and 



4 


/dX\ 







< 



1. (69) 



(d) ciJolAtl ~ 0, tuolAtl ^ 0. In this region, Eq.(|64D can never be a strictly vahd approximation 
because the right-hand side of Eq. (|6^ is zero at this point. Note however that the static 
approximation will be acceptable provided that the left-hand side of Eq. (|6^) is small at this 
point, which is true if Eq . (|68|) and Eq. (|69D are valid. 

Observe that Eqs. (|68|j69| ) are more stringent than Eq.(^). Evaluating Eq.(|68D at the 
resonance, we find 

(cos 26*0 — b + a) 



1^0 d 



2p 



< 



1. (70) 



For Eq.(|69D we only need to replace a — > —a in the above equation. Assuming that there is 
no accidental cancellation between the various independent terms, Eq. (|70D implies. 



,Sm^QbdT, < Tl ,da, < , Tl 2p dt , 

~ — ^ — ~ — (71) 

^ 2p T dt^ 2 ' 'rfT' ' 2 5m2dT'' ^ ' 



where we have used db/dT = 6b/T and recall that u^o = In deriving Eq.([7lD we have 

also neglected a term proportional to (cos2^^o — b + a) which is less stringent than Eq.QTI]) 
because (cos 26*0 — b + a) is just the resonance condition. The first condition in Eq . (|7l|) 
is satisfied provided that 

where we have set b = cos 26o (which leads to the most stringent condition) and we have used 
dT/dt ~ — 5.5T^/Mp. In order to evaluate the second condition in Eq. (|7lD , observe that 

da da dLy da , , 

dT dL,^ dT dT' ^ ^ 

Assuming that there is no accidental cancellation between the two terms on the right-hand- 
side of the above equation, the second term in Eq.(^) implies the following conditions at the 
resonance, 

aa, < ,r2^ 2p dt. , da dL,^. < T^^ 2p dt . 



'dT'^^ 2 5w?dT^' ^ dL^^ dT ^ ^ ' 2 6m^dT' 



(74) 
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Using da/dT ~ 4a/T, and a ~ 1, then the first equation above gives approximately the same 
condition as the first equation in Eq.(^). The second condition in Eq.(0) gives a condition 
on the rate of change of lepton number at the resonance. Expanding this equation out we 
find that 

1^1 ^ lH^ ' I . y'^^G'-^-^T^ ^ 4 X 10-^ ^ (75) 

' ' '2 dT2V2GFn^ 22^2 VMeVy MeV, ^ ' 

where we have used = 2C(3)T'^/7r^ ^ T^/4.1. Note that we have also used Eg . ([T0|) for the 
collision frequency. Thus, for example, if we are interested in studying the region where the 
lepton number is initially created, then a necessary condition for Eq. ( p3D to be approximately 
valid is that the resonance must occur for temperatures satisfying Eq. (|72D . From Eq.(P^) 
(with cos 26*0 — 1), this implies that 

W ~ 9 X 10^2 (5 X 10"^) eV^, (76) 

for Uf, — Vg {Vfji^T ~ ^s) oscillations. The creation of Ly^ must also satisfy Eq.(^) at the 
resonance. This condition should be checked when using Eq.(p3[) for self consistency. 

Perhaps surprisingly, there is a significant region of parameter space where the oscillations 
are not adiabatic at the resonance [i.e. 7 ~ 1 in Eqs.(|9],|0D] but Eqs.(O) 

are nevertheless 



satisfied. This is possible because Eqs. (|7lD are not equivalent to the adiabatic conditions 
Eqs.dliO). This is because Eqs.(|7l|) arise from demanding that the total contribution to 
dLy^/dt reduce approximately to the simple Eq. (|23|) . Recall that the total contribution to 
dLjj^/dt can be separated into two distinct contributions: from oscillations due to collisions 
and from oscillations between collisions. The adiabatic condition, on the other hand, is 
a necessary condition for the contribution of dL^^/dt from collisions to reduce to Eg . (|23|) . 
Thus it turns out that in the region when the system is both non-adiabatic and Egs.(|7T|) 
are satisfied, the modification to the equation for dLy^/dt from collisions which arises from 
the non-adiabaticity cancels with the extra contribution to dL^^/dt from oscillations between 
collisions. This type of cancellation is more transparent in the Hamiltonian formalism (see 
the appendix). 

Finally, to illustrate the analysis of this section, consider the examples given in Figures 1 
and 2. Recall that the solid and dashed lines correspond to the density matrix Eqs.(|^) and 
Eq.(^3]) respectively. Observe that for the example in Figure 1 (which has 5m?' = —1 eV'^), 
Eq.(^3]) is not a very good approximation at the resonance where the lepton number is initially 
created (although it is a reasonable approximation for small sin^ 26q) . This is because the 
lepton number is created so rapidly that Eq.(^) is not valid. However, for the example shown 
in Figure 2, where = —1000 eV"^, the temperature where the lepton number is created 
is much higher. Observe that Eq.(^) is not as stringent for high temperatures and it is 
therefore not surprising that the static approximation is approximately valid for this case. 
[Note that the result that the static approximation tends to be a good approximation at high 
temperatures can also be seen by observing that for high temperatures, ujq 0, and in this 
limit, Eq.(|64D will be satisfied]. 
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IV The Thermal Momentum distribution of the neutrino 



Hitherto we have made the assumption that the neutrinos are monochromatic. This as- 
sumption is not expected to hold for the neutrinos in the early Universe. The momentum 
distribution of these neutrinos will be the usual Fermi-Dirac distribution. Note that the width 
of the initial resonance in momentum space is much smaller than the spread of neutrino mo- 
menta. This means that only a few of the neutrinos will be at resonance at a given time. 
Also, not all of the neutrinos will be creating lepton number. Neutrinos in the region defined 
by > cos26'o (which includes part of the resonance) destroy lepton number, and those in 
the region < cos 26o create lepton number. The point where net lepton number is created 
only occurs when the lepton number creating neutrino oscillations dominate over the lepton 
number destroying oscillations. Recall that in the unphysical case where all of the neutrinos 
are assumed to be monochromatic, all of the neutrinos enter the resonance at the same time, 
where they all destroy lepton number if 6 > cos 26*0, or all create lepton number if 6 < cos26'o. 
Clearly, the effect of the thermal spread of momentum will make the creation of lepton number 
much smoother. An important consequence of this is that there will be even larger regions of 
parameter space where the system is smooth enough so that the static approximation is valid 
and hence Eq.(^) will be a good approximation (modified to incorporate the momentum 
dependence) . 

In the static limit, we can simply re-derive Eq.(^), assuming that the neutrino momenta 
form the usual Fermi-Dirac distribution. In this case, Eq.(^) becomes 



n. 



dL, 



dt —2 



\ J [T{iya Us) + r(z/„ Us)] {dn^^ - dn^^ - dn^^ + c/n^J, (77) 



2 



where 



_ 1 p^dp _ 1 p^dp 

~ 27r2 1 + e^P-^^)/T ' ~ 27r2 1 + e^P+^^)/^ ' ^ ' 

and driy^ , drtp^ are the differential number densities for the sterile and anti-sterile neutrinos 
respectively. In Eq.(^) /i is the chemical potential. 

The reaction rates can easily be obtained following a similar derivation as before [Eqns.(^ 
rSD], but this time we keep the momentum dependence (rather than setting p = (p)). Doing 



this, we find the following equation for the rate of change of lepton number in the static limit: 
dL.^ vr^ f s'Tiy{c-}f){dnt^-dnl) 

■ dt 4C(3)T3 J [xP + {c-bP + aP)2] [xP + {c-hP- aP)^] ' ^ ' 

where A is a small correction term 

-vr^ f s'TlJxP + {aP)^ + {hP-c)\dnZ^~dn-) 
8C(3)T3i [xP + (c - bP + aP)^][xP + [c - bP - aP)^] ' ^ ' 

and dn^^ = dn,,^ ± driy^. Recall that c = cos 26^0, s = sin26'o. In these equations note that 
the quantities, b^, qP, xP, Tp^ are all functions of momentum of the form: 

V P V 2 , ^ Ir. I P \ f '"^P V T^v -n P 



If = b-^ aP = a— a;P = + — r? = r J— 
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where a,b,T^^ are defined in Eqs.(|TH|, |TD|). Eq.([75|) reduces to Eq.(^), in the limit where all 
of the neutrino momenta are fixed to p = (p). [Note that aP = a when p = (p) ~ 3.15T and 
similarly for lP,x^ and T^^]. 

Note that the chemical potential is related to the lepton number by the equation 



Using Eqs.(|78|,|2|) we find 



dn} 



1 p^dp 
7r2 1+eP/' 



dn,. 



+ 0{LIJ. 



Thus substituting the above relations into Eq.(|75D, we find that 



dL, 



TT 



p" 



dn'! 



7r2(l + eP/r) dp 



dt 4C(3)T3 Jo [xP + {c-bP + aPy] [xp + {c - bP - aP)^ 

where A is a small correction term 

s^TP [xP + (aP)^ + {bP- c)2] ri2C(3)L^^p2gp/T 



A 



-TT 



8C(3)T3 Jo [xP + {c-bP + aP)2] [xp + {c - bP - aP) 



7r4(l + eP/^y 



dp 



(83) 



dp + A, 
(84) 



dp. 



(85) 

Eq.(p4D can be integrated numerically to obtain L^^ as a function of time (or temperature). 
We will give some examples in the next section after we discuss how to calculate the sterile 
neutrino number distributions. 

The main effect of the thermal spread of neutrino momenta is to make the generation of 
lepton number much smoother. From a computational point of view, this is very fortunate. 
This is because Eq.(^), like Eq. (P^D , is only valid provided the lepton number generation 
is sufficiently smooth (see the previous section for a detailed discussion of this point). In 
particular, Eq. (p^ should be a much better approximation to reality at the resonance where 
significant lepton number is initially generated. 

To complete this section, we comment on the rate of change of lepton number due to 
ordinary-ordinary neutrino oscillations. For definiteness consider Ue — oscillations. The 
rate of change of Ly^ — L^^ is given by 



dt - - I r('^/. ^ Jye)dny^ + J r(z/^ iye)dnp^ 
+ J T{ue i^ii)dny^ - J r(z/e h'fj)dnu^. 



(86) 



Using r(z/^ — > Ug) = r(z/e u^) (and similarly for the anti-neutrino rates), Eq.(p6D becomes 



di 



+ rr(z7. 



/o~r(z/^ 

^e)( 



V, 



driuf^ 
dp 



dui, 



dp 



dp 



dp 



dp 



dp, 



11) 



20 



where 



dn^ 1 p2 dna 1 



dp 27r2 l+e{P-A'i)/T ; dp i+e(P+Mi)/T ) 

dp 27r2 l+e(P-A'2)/r' dp 27r2 l+e(P+M2)/T • 

The chemical potentials /ii 2 are related to the lepton numbers L^^ ^ by the equations 

where we have assumed that fii/T <^ 1. Using these relations and expanding out Eq. (|8^) 
(again assuming that fii/T <^ 1) we find to leading order that 

~ ' ' [r(z/^ Z/e) + r(z/^ ^ Z/e)] c/p. (90) 



dt 7r2T3 Jo (1 + eP/^) 

From the above equation we see that L^^ — L^,^ always evolves such that [L^^ — L^J 0. 
Also note that Eq.(PD|) shows that the rate of change of lepton number due to ordinary- 
ordinary neutrino oscillations is generally smaller than the rate of change of lepton number 
due to ordinary - sterile neutrino oscillations (assuming L^,^ <^ 1). [Actually Eq.(^) has a 
strength comparable to the correction term A for ordinary-sterile neutrino oscillations Eq. (p^) , 
although note that the mixing angle between ordinary neutrinos can be significantly larger 
than the mixing angle between ordinary and sterile neutrinos]. 

For ordinary-ordinary neutrino oscillations, neutral current interactions do not collapse the 
wavefunction because they cannot distinguish different neutrino species. Only the charged 
current interactions can distinguish the neutrino species. For example, for temperatures 
1 MeV ~ T ~ 30 MeV, there are near equilibrium number densities of electrons and positrons. 
The number of muons and anti-muons will be much less than the number of electrons and 
positrons, and we will neglect them (actually these are important for z/,- — z/^ oscillations). The 
rate at which charged current interactions occur is given approximately by Fj, ~ iFj,^ ~r,y^| ~ 
y'^G^pT^, where y'^ ^ i/e — [— 1.1 see Eq.(0)]. Also note that anti-neutrino -neutrino pO| 
and neutrino - neutrino ||31|| forward scattering amplitudes induce off diagonal contributions to 
the effective potential. [Note that these contributions do not occur for the effective potential 
governing ordinary - sterile neutrino oscillations]. It would be necessary to include these 
effects in-order to evaluate the reaction rates. We leave this as a take home exercise for the 
reader. 



V. The effects of non-negligible sterile neutrino number densities and the param- 
eter space for large lepton number asymmetry generation 

In this section we will do three things. We will study the effects of non- negligible sterile 
neutrino number densities, which can arise for the case of relatively large, or even moderate 
values of sin^ 26*0. We will examine the parameter space where significant generation of lepton 
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number occurs. Finally, we will obtain the BBN bound on the parameter space for ~ 
oscillations with Sm? < 0, and with |5m^| ~ 10"'' eV^, sin^ 26^0 ~ 10"^. 

There are several ways in which the creation of lepton number(s) can prevent the sterile 
neutrinos from coming into equilibrium. One way is that one set of oscillations z/^ — creates 
L^^ . The lepton number L^^ can then suppress other, independent oscillations such as — Vs 
oscillations (with (3 ^ a) for example. A more direct, but less dramatic way in which the 
creation of lepton number can help prevent the sterile neutrinos from coming into equilibrium, 
is that the lepton number generated from say — oscillations itself suppresses the — 
oscillations ||32||. We will examine the latter effect here (some examples of the former effect 



will be studied in the next section). Previous work]^, ^, 0, |^ obtained the BBN bound for 
large \5m?\ ~ 10~^ eV'^ (with 5m? < 0) and small sin^ 26o ~ 10^^ which can be approximately 
parametrized as follows p[ 

|W|sin^2^o ~ 10"^ e1/l (91) 

This bound arises by assuming that the — Vg oscillations do not bring the sterile Vg state 
into equilibrium. Note that this bound neglected the creation of lepton number and it also did 
not include the effects of the distribution of neutrino momentum. However, in the realistic 
case, the creation of Ly^ (after it occurs) will suppress the Va — J^s oscillations and the actual 
bound would be expected to be somewhat less stringent than Eg . (pTj) . 

To proceed we will need to examine the effects of non-negligible sterile neutrino number 
densities. The evolution of the number distribution of sterile neutrinos is governed by the 
rate equation 

rl rln /rln rin / rin — rlrt / rin 

I/, ^ Us). (92) 
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driyjdp 




dn^^/dp - dn^Jdp 


dt 


_dny^/dp_ 




driy^/dp 



A similar equation holds for the number distribution of sterile anti-neutrinos. Introducing 

djii/j /dp \ 
drii)^ /dp)'' 



the notation, z = '^""^ ^ff (for anti- neutrinos we use the corresponding notation, z = *liisMi' 



Eq.(p2D becomes 



— = {\-z)Y{v^-^Vg) = ^ ^ -. 93 



The corresponding equation for anti-neutrinos can be obtained by replacing z z and 
—aP in the above equation. In solving the above differential equation, we will assume the initial 
condition ^ = 0. We will also assume that the number densities of the ordinary neutrinos 
are given by their equilibrium values. Note that the quantity z depends only on the reaction 
rates and is otherwise independent of the expansion. 

From the definition of z, it follows that dug = zdn,y^, dug = zdrip^. Thus, from Eq. 



dLy. ^. 1 r s^Tiy{c-}f) {l-z+)p'dp ^g^^ 



dt AC{3)T^Jo [xP + {c - bP + aPy][xP + {c ~ bP - apy] 1 + eP/^ 
where A is a small correction term 

1 s'rPjxP + {aP)' + {bP-c)'] z-p'dp 

8C(3)r3io [xP + {c-bP + aPy][xP + {c-bP-aPy]l + eP/^' ^ ' 
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with = {z ± z)/2 and we have neglected a small correction term which is proportional to 
L^^. Note that Eq.(p3D and Eq.(0) must be solved simultaneously. 

For the numerical work, the continuous variable p/T is replaced by a finite set of momenta 
Pn/T (with n = 0,1,. ..,A^) and the integral over momentum in Eq.(p^) is replaced by the 
sum of a finite number of terms. Correspondingly, the variable z{t,p/T) is replaced by the set 
of variables, Zn{t), where the evolution of each variable Zn{t) is governed by the differential 
equation, Eq . (|93D [with p/T = Pn/T for z = z^it), n = 0, 1, N]. Thus, the single differential 
equation, Eq.(p3D is replaced by a set of differential equations, one for each momentum 
step. These differential equations, together with Eq.fplf), are coupled differential equations 
which must be integrated simultaneously. 

We now illustrate the creation of lepton number as governed by Eqs.(0,0) with some 
examples. We have numerically integrated Eqs. ( P^^BD for the following parameter choices. In 
Figure 3 we have considered u^^t — Vg oscillations with the parameter choice 5m^ = — 1 eV^^, 
sin^ 2^0 = 10"^ (dashed line), sin^ 26^ = IQ-^ (dashed-dotted line) and sin^ 2^o = 10"*^ (solid 
line). Figure 4 is the same as Figure 3, except that Sm^ = —1000 eV"^, sin^ 29o = 10^^ (dashed 
line), sin^ 2^0 = 10~^ (dashed-dotted line) and sin^ 26'o = 10~^ (solid line). In both examples 
we have assumed that the initial lepton asymmetry is zero. Recall that the generation of 
lepton number is essentially independent of the initial lepton number asymmetry provided 
that it is less than about 10~^ (for more discussion about this point see section 2). Note that 
for convenience we have plotted \Ly^\. The lepton asymmetry L^^ changes sign at the point 
where it is created. Before this point L,^^ evolves such that it has the opposite sign to t] while 
for evolution subsequent to the point where L^^ is initially created, L^,^ has the same sign as 
Tj. Recall that this behaviour is expected (see the earlier discussion in section II). 

In these examples, the generation of lepton number is considerably smoother than in the 
earlier case where the momentum distribution was neglected (see Figures 1,2). For this reason, 
it turns out that throughout most of the evolution of L,y^ , the rate of change of L^^ satisfies 
the condition Eq.(^) and thus Eq . (|8^) should be approximately valid (except at quite low 
temperatures where the MSW effect will be important). 

In order to gain insight into the effects of the neutrino momentum distribution, it is useful 
to compare Figures 3,4 (which incorporate the neutrino momentum distribution) with the 
Figures 1,2 (where all of the momentum of all of the neutrinos were set equal to the mean 
momentum). Qualitatively, there is not a great deal of difference. However there are several 
very important effects, which we summarize below. 

(1) For the examples with relatively small sin^2^o, lepton number creation generally begins 
somewhat earlier (i.e. at a higher temperature) than in the case where momentum distribution 
is neglected. For the examples shown in Figures 3,4 with 6m'^ = —1 eV"^, sin^ 2^o = 10~* 
(-1000 e1/2, sin2 2^o = 10"^), lepton number is created when T ~ 20 MeV (T ~ 65 MeV). 
This can be compared with the simplistic case where the neutrino momentum distribution 
was neglected. In this case, we see from Figures 1, 2 that lepton number creation begins at 
T ~ 16.0 MeV (T ~ 50 MeV) for 6m^ = -1 eV^ (-1000 eV^). The fact that the critical 
temperature increases can be explained rather simply. Note that the neutrino number density 
distribution peaks at about p ~ 2.2T, which should be compared with the average momentum 
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of about 3.15T used in Figures 1 and 2. Using the former approximation instead of the latter 
leads the critical temperature to increase by about 12%. This explains qualitatively why 
the critical temperature increases. Note, however, that the accurate numerical calculations 
displayed in Figures 3 and 4 actually show that the temperature increases by more than this, 
and also that the temperature increase depends on the mixing angle. 

(2) For the examples with large sin^ 29o, the point where significant generation of lepton 
number is created occurs much later than in the examples with small sin^ 29q. The reason for 
this is that for large sin^ 29q, the number density of sterile neutrinos is larger. In the region 
before significant lepton number is generated, a ~ and all of the neutrino oscillations with 

< cos 2^0 have already passed through the resonance while the neutrino oscillations with 
> cos26'o have yet to pass through the resonance. Since the creation of sterile neutrinos 
is dominated by the oscillations at the resonance, it follows that the sterile neutrino number 
distribution with momenta in the region where b^ < cos2^o will be much greater than for 
sterile neutrinos with momenta in the region where b^ > cos2^o- Thus, from Eg. (p4D , the 
lepton number creating oscillations (with b^ < cos2^o) are suppressed if the number density 
of sterile neutrnos is non- negligible, as occurs for large sin^ 29q. The lepton number destroying 
oscillations (with b^ > cos2^o); on the other hand, are not suppressed because the number of 
sterile neutrinos with b^ > cos 29q are negligible. 

(3) The creation of lepton number is considerably smoother in the realistic case. For instance, 
in the example where = —1 eV"^, sin^ 26'o = 10~^, in the realistic case (shown in Figure 
3), Lj,^ goes from 10~^° to 10~^ in about AT ^ 1 MeV, whereas in the unrealistic case where 
the neutrino momentum distribution was neglected (shown in Figure 1), \L^J goes from 10~^° 
to 10-6 in about AT ^ 0.005 MeV. 

(4) At low temperatures, the lepton number gets "frozen" at an earlier time. For example, in 
the case where Srri^ = —1 eV"^ and sin^ 26'o = 10"^, with momentum dependence (Figure 3), 
the final value for the lepton number is ~ 4 x 10"^, whereas in the unrealistic case without 
the neutrino momentum distribution, the final value for the lepton number for this example 
(Figure 1) is ~ 10~^. As discussed briefly in section II, this effect is expected because the 
temperature where the lepton number gets frozen occurs when the rate of change of the 
variable a due to the expansion of the Universe dominates over the rate of change of a due 
to neutrino oscillations. In the realistic case where the momentum distribution is taken 
into account, the maximum value of the rate of change of a due to neutrino oscillations is 
suppressed because only a small fraction of the neutrinos will be at the resonance. 

This last point suggests that the momentum distribution cannot be ignored if one is 
interested in finding out the precise final value of the lepton number generated. However, note 
that Eq.(p4D does not incorporate flavour conversion due to the MSW effect [see assumption (4) 
in section II for some discussion about this point]. The effect of the MSW flavour conversion 
should be to be to keep a ~ 1 for lower temperatures. This means that the final value of 
Ly^ should be significantly larger than suggested by Figures 3,4. This effect will need to be 
incorporated if one wants to calculate the precise value of the final lepton number generated. 
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[The precise value of the final lepton number can be obtained by numerically integrating the 
density matrix equations Eqs. (^61) suitably modified to incorporate the neutrino momentum 
distribution]. In particular, if one is interested in working out the region of parameter space 
where the electron lepton number is large enough to affect BBN through nuclear reaction 
rates, then the final value of the electron lepton number is very important |]33|, |34| . 

Note that we can check Eq.(|9T|) by numerically integrating Eq. (p3|) and Eq.(p4D assuming 
for definiteness that PuJ Pu^ ~ 0.6 (where the p's are the energy densities). This leads to the 
following constraint on 5m^, sin^ 26'o, 



sin2 2^o ~ 2(4) X 10" 



(96) 



for Ve — i^s {i^ii,T — Vs) oscillations. Thus, we see that Eq.(|^) turns out to be a good approxima- 
tion after all. This is basically due to the result that the creation of a non-negligible number 
of sterile neutrinos has the effect of delaying the point where significant lepton number is 
created [see point (2) above]. 

Finally, the region of parameter space where significant neutrino asymmetries are gen- 
erated by ordinary sterile neutrino oscillations can be obtained by integrating Eq. (p4D and 
Eq. (P^D . The result of this numerical work is that significant neutrino asymmetry (|-^vi^^| ~ 
10~^) is generated by ordinary - sterile neutrino oscillations for the following region of param- 
eter space 



6(5) X 10 



-10 



\5m?\ 



~ sin2 2^o ~ 2(4) X 10" 



\5711? 



and \5rn? 



~ 10" 



(97) 



Ug) oscillations. Note that we have assumed that PuJ Pua ~ 0.6 [Eq. (p6| 



for z/e 

In the general case where no bound on PuslPvc is assumed, the upper bound on sin^ 26*0 is 



considerably weaker. For example. 



oscillations with 5m 



■1 eV' 



sm' 



2^n = 10-^ 



violate the bound, Eq.(p6D but still generate significant neutrino asymmetry, as illustrated in 
Figure 3. [For this particular example, we found that p^Jpu^ — 0.86]. 

The parameter space in Eq. ( pTf ) can be compared with previous work where the momentum 
dependence was neglected jlSl, [l2| . As we have mentioned above, the upper bound on sin^26'o 
which assumes a BBN bound of PvJ Pv^ ~ 0.6, is not modified much when the momentum 
distribution of the neutrino is incorporated. For the lower limit of sin^2^0; the effect of the 
momentum dependence is to reduce the region of parameter space by nearly two orders of 
magnitude. 

Finally, it may be possible for significant neutrino asymmetries to be generated for |5m^ 
10~^ eV"^, however the mechanism of production of these asymmetries is dominated by 
oscillations between collisions (rather than the mechanism of collisions) and tend to be 
oscillatory [l^ 



< 
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VI. Consistency of the maximal vacuum oscillation solutions of the solar and 
atmospheric neutrino problems with BBN 
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We now turn to another application of the phenomenon of lepton number creation due to 
ordinary-sterile neutrino oscillations. First, in the context of a simple explanation of the 
solar neutrino problem which involves large angle Vf, — Ug oscillations, we will determine the 
conditions under which the lepton number produced from Ua — i^s oscillations can suppress the 
oscillations — Ug (where (3 a). This allows the BBN bounds on ordinary-sterile neutrino 
oscillations to be evaded by many orders of magnitude, as we will show. We begin by briefly 
reviewing the maximal vacuum oscillation solution to the solar neutrino problem 0. 

One possible explanation of the solar neutrino problem is that the electron neutrino os- 
cillates maximally (or near maximally) with a sterile neutrino (which we here denote as u'^ 
rather than as Ug in order to remind the reader that this sterile neutrino is approximately 
maximally mixed with t'e)[0, |3^- We will denote the Sm^ for z/g — z/g oscillations by Sml^,. As 
is well known, for a large range of parameters 

3xlO~^°eV^2~ I^^L'I- 10"^ e1/^ (98) 

maximal vacuum oscillations imply that the flux of electron neutrinos from the sun will be 
reduced by a factor of two for all neutrino energies relevant to the solar neutrino experiments. 
We will call this scenario the "maximal vacuum oscillation solution" to the solar neutrino 
problem. It is a very simple and predictive scheme which can either be ruled out or tested 
more stringently with the existing experiments. Importantly, it also makes definite predictions 
for the new experiments, SNO, Superkamiokande and Borexino. Our interest in this scheme is 
also motivated by the exact parity symmetric model (see Ref.|]T^ for a review of this model). 
This model predicts that ordinary neutrinos will be maximally mixed with mirror neutrinos 
(which are approximately sterile as far as ordinary matter is concerned) if neutrinos have 
mass[]T^. If we make the assumption that the mixing between the generations is small (as 
it is in the quark sector) then the exact parity symmetric model predicts that the three 
known neutrinos will each be (to a good approximation) maximal mixtures of two mass 
eigenstates. There are also other interesting models which predict that the electron neutrino 
is approximately maximally mixed with a sterile neutrino p7||. The maximal mixing of the 
electron neutrino (z/g) and the sterile neutrino will reduce the solar neutrino flux by an energy 
independent factor of two for the large range of parameters given in Eg . (p8|) . This leads to 
definite predictions for the expected solar neutrino fluxes for the existing experiments. In 
Ref. P], we compared these predictions with the existing experiments. We summarize the 



results of that exercise in Table 1 which we have updated to include the most recent data [p^] . 

Note that in Table 1, the Kamiokande experiment has been used as a measurement of the 
Boron flux[p9|, |40| , p| . This is a sensible way to analyse the data (but not the only way of 



course) because the flux of neutrinos coming from this reaction chain is difficult to reliably 
calculate P2|. Clearly, the simple energy independent flux reduction by a factor of two leads 
to predictions which are in quite reasonable agreement with the data. If the minimal standard 
model had given such good predictions, few would have claimed that there is a solar neutrino 
problem. 

Note that the maximal vacuum oscillation solution is distinct from the "just so" large angle 



vacuum oscillation solution|H5|. In the "just so" solution, the electron neutrino oscillation 
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length is assumed to be about equal to the distance between the earth and the sun (which 
corresponds to \5m?\ ~ 10~^° eV"^). In this case the flux of neutrinos depends sensitively on 
Sw? and it is possible to fit the data to the free parameters dm?, sin^ 2^o|0- The advantage of 



doing this is that a good fit to the data can be obtained (however this is not so surprising since 
there are two free parameters to adjust). The disadvantage is that fine tuning is required and 
predictivity is lost because of the two free parameters. The maximal mixing solution on the 
other hand assumes maximal mixing and that is in the range Eg . (pSf ) . For this parameter 
range there is an energy independent flux reduction by a factor of two. The advantage of this 
possibility is that it does not require fine tuning and it is predictive. A consequence of this 
is that it is testable with the existing experiments. The disadvantage of this scenario is that 
it does not give a perfect fit to the data. However, in our opinion the predictions are in 
remarkably good agreement with the data given the simplicity and predictivity of the model. 



With the range of parameters in Eq.(^) there is a potential conflict with BBNP3|, ^ 
For maximally mixed z/g and z/g neutrinos, the following rather stringent BBN bound has been 
obtained assuming that the lepton number asymmetry could be neglected^, |^, |^, |[ : 

\6ml^,\ ^ 10''' eV\ (99) 

This bound arises by requiring that the sterile neutrinos do not significantly modify the 
successful BBN calculations. For temperatures above the kinetic decoupling temperature 
the requirement that the sterile neutrinos do not come into equilibrium implies the bound 
|5mgg,| ~ 10~^ eV^. Smaller values of Sml^, in the range 10~® ~ \6ml^,\/eV^ ~ 10"^ can be 
excluded because the oscillations deplete the number of electron neutrinos (and anti-neutrinos) 
after kinetic decoupling (so that they cannot be replenished). The depletion of electron 
neutrinos increases the He/H primordial abundance ratio. This is because the temperature 
where the ratio of neutrons to protons freezes out is increased if there are less electron neutrinos 
around. For ~ 10^*^ eV"^, the oscillation lengths are too long to have any significant 

effect on the number densities of electron neutrinos during the nucleosynthesis era. If the 
bound in Eq. (|99D were valid then it would restrict much of the parameter space for the maximal 
vacuum oscillation solution of the solar neutrino problem. However, this bound does not hold 
if there is an appreciable lepton number asymmetry in the early Universe for temperatures 
between 1 — 100 MeV[]TT|. This is because the generation of significant lepton number L*^^^ 



implies that the quantity age' [which is the a parameter defined in Eq.(|18D with = Sml^,^ 
is very large thereby suppressing the oscillations [note that for a^e' ^ 1, sin < sin^ 2^0, 
see Eq.(0)]. We will now show in detail how the creation of lepton number can relax the 
BBN bound Eg . (|99D by many orders of magnitude. 

We will assume that the various oscillations can be approximately broken up into the 
pairwise oscillations z/g — z/g, I'fj, — z/g and Ur — z/g. We will denote the various oscillation 
parameters in a self-evident notation, 

bae',o,ae' for z/q, — Z/g oscillatious, (100) 

where a = e,fi,T. We will denote the mixing parameters, (5m^, sin^ 2^o appropriate for 
fa — f'e oscillations by Sm"^^, , sin^ 29^'^ . Note that lepton number cannot be created by 



27 



Ua — i^'e oscillations until b^e' < cos26'q^'. Recall that the b parameter is inversely proportional 
to [see Eq. ([l8|) ]. Thus, the earliest point during the evolution of the Universe where lepton 
number can be created due to ordinary - sterile neutrino oscillations occurs for oscillations 
which have the largest |5m^|. Note that these oscillations must satisfy the bound in Eq.(|9 
if they are to produce lepton number, and they should also satisfy the BBN bound Eq.(|9 
if we demand that the sterile neutrino energy density be small enough so that BBN is not 
significantly modified. Note that the J^e — K oscillations have very small|5m2^,| ~ IQ-^ eV^§^, 



and cos (assuming maximal or near maximal mixing), and thus these oscillations 

themselves cannot produce significant lepton number. However, the 6m'^ for Ur — z/g or — v'^ 
oscillations can have much larger 5m^ (and they should also have 5m? < if niy^^my^ > 
my/)[0. We will assume for definiteness that m^^ > m^^ > mu'^ so that |5m^g, | > l^m^e'l 
and the z/,- — u'^ oscillations create Ly_^ first (with Ly^,Ly^ assumed to be initially negligible). 
If rrii,^ > TTLy^ then we only need to replace z/^ by z/^ in the following analysis. 

Thus, we will consider the system comprising z^r, J^e and z/g (and their anti-particles). 
Our analysis will be divided into two parts. First, we will calculate the condition that the 
L^'^'> created by Ur — v'^ oscillations survives without being subsequently destroyed by z/g — z^e 
oscillations. We will then establish the conditions under which L*^^^ is created early enough 
and is large enough to suppress the z/g — v'^ oscillations so that only a negligible number of v'^ 
is produced. 

For simplicity we will first analyse the system neglecting the momentum distribution of 
the neutrino. This is useful because under this assumption it turns out that this system can 
be approximately solved analytically as we will show. We will then consider the realistic case 
where the spread of momenta is taken into consideration. 

It is important to observe that the generation of Ly^ also leads to the generation of 
L(^) [through Eq.(|l6D]. If we assume that negligible Ly^ is generated, then L^^^ ^ L^'^'^ /2. 
However z/g — z^g oscillations can potentially generate Ly^ such that L^^^ 0. (Recall that 
is an approximately stable fixed point for the z/g — z/g system for temperatures 
greater than a few MeV). The effect of the z/g — z^g oscillations will be greatest when the 
z/g — z/g oscillations are at resonance. If negligible Ly^ is generated, then |age'| ~ R\aT-f.i\/2 
and |6gg'| ~ R{A(,/ Ar)\bT-ei\ (where R = \6m'^^, / 6ml^,\) . Hence the z/g — z/g resonance condition 
(ogg/ = 6gg'), will be satisfied when 

\are'\=2{AJA^)\bre'\. (101) 

Recall that the Ur — z/g oscillations generate Ly^ such that 

are' ^ 1 - bre', (102) 

where we have assumed that cos 2^™' ~ 1 and that Ly^ > for definiteness. Observe that 
Eqs.( |101] , |102|) imply that the system inevitably passes through the z/g — z/g resonance. This 



event will occur when 

|6,g,|~- (103) 

' ' + 2Ag ^ ^ 
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Using the definition of bre' which can be obtained from Eq.(|TH|), the above equation can be 
solved for the z/g — i^'^ resonance temperature 



res 




MeV. (104) 



Thus, when T = T^^'si ^^e z/^ — z/^ oscillations have created enough L*^*^) so that the — ^'e 
oscillations will be at the resonance, assuming that negligible L^^ has been generated. In 
general the z/e — z/g resonance temperature depends on both L^^ and Lj,^ . The z/g — z/g resonance 
condition Oee' = &ee' implies that the resonance temperature for z/g — z/g oscillations is related 
to Lj^^ and Ly^ by the equation 



T:i = ^'^m^^^{2L^^ + Ul (105) 

where we have neglected the small baryon and electron asymmetries and a possible mu neu- 
trino asymmetry (we will discuss the effects of the mu neutrino later). Thus the resonance 
temperature will change when L^^ and Ly^ change due to oscillations. 
Let us consider the rate of change of the quantity {T^^^ — T), 

d{T^i - T) _ gT-: dU^ dT^i dLy^ dT 

dt dLy^ dt dLy^ dt dt' ^ ^ 

evaluated at the temperature T = T^e's- Note that the first term on the right-hand side 
of Eq. (|106|) represents the rate of change of T^^^ due to Ue — J^'e oscillations (and its sign is 
negative), while the second term is the rate of change of T^^^ due to Ur — J^'e oscillations (and 
the sign of this term is positive). The third term in Eq.( |106|) is the rate of change of [T^gg ~T) 
due to the expansion of the Universe (-^ — 5.5T^/Mp) (this term is also positive in sign). 
Observe that if d{T^^'^ — T)/dt > 0, then the system passes through the resonance without 
significant destruction of L^^\ If on the other hand, d(T^^^ —T)/dt ~ 0, then the position 
of the resonance moves to lower and lower temperatures and L*^"^) — >■ 0. Thus, a sufficient 
condition that L^^^^ survives without being destroyed by Ve — v'^ oscillations is that 

dL,^ dt dL,^ dt dt' ^ ' 

To evaluate dL^^/dt, observe that 

dLy dLy -AL^^'^ dT , , 

where we have assumed that L^'^^ ~ T~^ for T = T^l^. Of course this latter assumption only 
holds provided that the z/g— z/g resonance does not occur while Ly^ is still growing exponentially. 
However, for sin^ 2Qq^ sufficiently large, the z/g — z/g resonance can occur during the rapid 
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exponential growth phase of L^^ . If this happens then the rate at which the Vj. — v'^ oscillations 
move the system away from the — z/^ resonance is much more rapid. Consequently, the 
region of parameter space where L^^^ survives without being destroyed by z/g — v'^ oscillations 
is significantly larger in this case (this effect will be illustrated later on when we study the 
system numerically). 

Thus using Eq.(|108|), Eq. (|107|) can be written in the form 



3 5L„„ dT: 



4 dt dU 



res > 



1 + 



res 

T 9L,„ 



dT 



(109) 



where we have used the relation dT^^'^/ dLy^ = 2dT^^'^/ dL^^ which is easily obtainable from 
Eq.(|Ol). 

Note that the most stringent condition occurs at the z/g — z/g resonance temperature, 
Eq.(104). We are primarily interested in relatively large values of |5wi^e'l ~ 10"^ eV'^, which 
means that T^^'^ ~ 8 MeV. Thus, from section III, we are in a region of parameter space 
where we expect Eq. (p3D to be vahd. [In particular, note that since T^^^ is not at the point 
where the lepton number is initially created, Eq . (|75|) should also be valid]. Thus, from Eq. (|2^) 
we can obtain the rate of change of L^,^ due to Ue — v'^ oscillations, at the z/g — v'^ resonance 
(where 6 — a — c = 0). We find 



dt 



3 sin^ 9f 



6.3 



(110) 



Note that from Eq. (|105|), we have: 



res 



Ml 

A ''Pee' 
^e-'- res 



r-pe 



'III) 



Thus, 



Tee' f)T 
-^res '-'■^i^e 



(112) 



Hence the sufficient condition that L^^'^ survives without being destroyed by V(. — v'^ oscillations 
can be obtained by substituting Eqs. (|110| , |lll| , p.l2| ) into Eq. (|109|) . Doing this exercise we find 



1 5m; 



~ A|(5m: 



2 I — 

re' I ) 



(113) 



where A is given by 



A 



12.6Gf 



Y^sin^ 2Qlf Mw 
12.6Gf 



/ 8?/eAe5.5 \ 

V 3Afp } 



I (|5m2|)i/6 



/ 8i/e^e5.5 \ 

V 3Mp j 



4.1M2, 



(114) 
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and we have used Eq.( |104| ). Thus, putting the numbers in, we find 

(115) 

where we have assumed maximal mixing (i.e. sin^ 26'q'^' ^ 1). Thus provided this condition 
holds, L^^^ will not be destroyed significantly by Ue — v'^ oscillations (under the assumption 
that the neutrino thermal momentum distribution can be neglected; we will study the effects 
of the momentum distribution in a moment) and the system moves quickly away from the 
Vf. — v'^ resonance. While this condition was derived as a sufficient condition, it turns out to be 
a necessary one as well. This is because if Eq. (|107| ) where not valid, then v^. — v'^ oscillations 
would create Ly^ rapidly enough such that diT^l^ - T)/dt < 0. This would mean that the 
Ve — ^'e resonance temperature would move to lower and lower temperatures where the rate 
of change of Ly^ [from Eq.( |llCI| )] would be even larger (as it is proportional to 1/T^) and the 



expansion rate slower. Thus if the condition Eq.( |107| ) were not satisfied initially, it could 
certainly not be satisfied for lower temperatures. 

In the above system consisting of z/,-, z/g and v'^ (and their anti-particles) discussed above, 
observe that we have neglected the effects of — oscillations. As discussed in the previous 
section, the effect of these oscillations is to make (L^^ — L^^) tend to zero. Since these 
oscillations cannot prevent L^^ from being generated the effect of incorporating them should 
only increase the allowed region of parameter space [^. 

The effect of the muon neutrino can also only increase the allowed region of parameter 
space. The effect of i/^ — z/g oscillations will be to create Ly^ provided that 5m^g, < 0. The 
effects of Pfj, are completely analogous to the effects of the z/^ neutrino, and we can replace Pr 
with z/^ in the above analysis. This means that it is only necessary that either or 6m?^^, 

(or both) satisfy Eq.( |115| ). 

Hitherto we have examined the system neglecting the thermal distribution of neutrino 
momenta. We now study the realistic case where the thermal distribution of neutrino momenta 
is taken into consideration. We first estimate approximately the effects of the momentum 
distribution analytically and then we will preform a more accurate numerical study. 

The previous calculation assumes that all of the neutrinos have a common momentum and 
thus they all enter the resonance at the same time. In the realistic case, only a small fraction 
(less than about 1 percent as we will show) of the neutrinos are at resonance at any given 
time. Note that the v^- — v'^ oscillations are not affected greatly by this consideration, since as 
we showed in section V, the momentum spread does not prevent Ly^ from being created (and 
it still satisfies approximately L'-^-' ~ after it is initially created). On the other hand the 
effect of the neutrino momentum distribution on the z/g — v'^ oscillations is very important. 
This is because the z/g — v'^ oscillations cannot destroy L*^^^ as efficiently as before. In fact 
Eq.(rTr) will be reduced by a factor which is about equal to the fraction of neutrinos at 



the resonance. In principle, one should solve Eq.( |109| ) at the point where electron neutrinos 
of momentum p = 7T are at resonance and then calculate the minimum of the value of A 
[see Eqs.( [L13"| , |114| )] over the range of all possible values of 7. For simplicity we will make a 



rough approximation and assume that the minimum of A occcurs when neutrinos of average 



31 



momentum are at resonance (i.e. assume 7 ~ 3.15). [Note that later-on we will do a more 
accurate numerical calculation]. 

To calculate the fraction of neutrinos at the z/g — v'^ resonance we need to calculate the 
width of the resonance in momentum space. We will denote this width by Ap. From Eq.p^) 
it is easy to see that the width of the resonance is governed approximately by the equation 



dp 



2v^ 



ee' ; 



(116) 



where we have assumed maximal mixing (i.e. cos 2^0*^' ~ 0). Note that from the momentum 
dependence of V^a^ [see Eq. (Pl)], it follows that 



and hence 



dp p ^ dp 



p 



P 



p 



(117) 



:ii8) 



dp p 

where we have used the result that 6gg, ~ a^g, at the resonance (note that we have assumed that 
L^^'^ > for definiteness) . Note that we are essentially interested in evaluating the maximum 
value of the fraction of neutrinos at the resonance. This maximum fraction should occur 
approximately when p ~ (p). Thus, from the previous analysis, the resonance for neutrinos 
of average momentum occurs when 



are' Sml^, ^ 1 5ml^, 
2 (5mgg/ 2 5m^g, ' 



(119) 



Thus, from Eqs.( |ll^ , |119| ) and Eq.( |116| ) the width of the resonance in momentum space 
becomes 

\5ml^,\ 



Ap ^ ApJiXee') 



Recall that {xee') is defined in Eq. (pSj) , and is given by 

(. ,) - sin^ 29-' + - sin^ 29^' + ^-^'^(3-^^)^^" 

{Xee') - sm 29, + 1 _ sm 2U, + ^^^^^^^2 • 

Expanding out {xee') at the temperature T ^ T^f^ [defined in Eq. ( |104| )], we find 



(120) 



(121) 



{Xee')=Sm^29'/ + 



yeGpM^A.l Ar 
2^f2Ar Ar + 2Ae 
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5mL, 
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ee' 



(122) 



where the last part follows provided that \Sm^^,\ ~ 10 ISm^^,]. Thus, using the above 
equation, Eq. (|120D simplifies to 

Ap 2yeGpM^A.l Ar 



p 



y/2Ar Ar + A^ 



1.4 X 10" 



(123) 
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Thus it is clear that only a small fraction of neutrinos will be at the resonance. We denote 
the fraction of electron neutrinos at the z/g — v'^ resonance by Any/riy. Note that IS^n^jny is 
given approximately by the equation 



Using Uy = f C(3)TV7r2, and 



we find 



Til, Tiy dp 
driu 1 



dp 27r2l + eP/^' 



(125) 



An, ^ 2y,GpM^^A.l {p/Tf ^ 1.4x10-^ {p/Tf ^ 

''"'^^'^ v^A, A, + 2Ael.5C(3)(l + e<f)/^) - 1.5C(3) 1 + e(f)/^ " ' 

(126) 

The effect of the momentum spread is thus to reduce the number of neutrinos at the resonance 
by the above factor. Multiplying Eq.( [LlCI| ) by this fraction and repeating the same steps which 

lead to Eq. (|115| ) we find that Eq. (|115| ) is weakened by the factor ^riy/riu ~ 10^""^. In other 
words the effect of the neutrino momentum distribution is to increase the allowed region of 
parameter space for which — ^'e oscillations do not destroy the L^^'' asymmetry created by 
Vj. — v'^ oscillations. This region of parameter space is given approximately by 

\S^le'\^^, (127) 



where A is given in Eq. (|114|) and JAi^/riu is given in Eq. (|126|) . Putting the numbers in, the 



above condition can be written in the form 



I'^"^-I^6xl0-('M^V^ (128) 



e1/2 \^ e\/2 

We now check this result by doing a more accurate numerical study of this problem. 

The rate of change of L^,^ and L^^ due to the Ur — ly'^, z/g — z/g oscillations can be obtained 
from Eq. (|9^) . This leads to the following coupled differential equations 

dL^ _ 1 TOO sin2 2eg^'rg^a;:^,(cos2eg^'-fe^^ ) {l-z+)p^dp 

dt - 4C(3)T3 Jo [^P^^+^^os2ef-lf^,+al,mxlMcos29f~bl~al,)^] (1+eWT) 

+ 8C{3)T3J0 ,+(cos2ef -b" ,+a'' ,)2][x^ ,+{cos2e9/ -b" ,)^] 1+eP/T , 

dL^ ^ 1 roo sin^ 261-' Tl^al^, (cos 291^' -IF^,) (l-z+)p^dp 

— 4C(3)T3 JO UP +(cos2e-<^'-6P ,+0^ ,)2][a;P ,+(cos 2611=' -6'' ,-0^ ,)2i (1+eP/T) 

L ^g/ \ U ^g/ ^g/ / J L \ y / J 

1 00 sin2 2e5'='rg^[x;^,+(a;^,)2+(b;^,-cos2e^'=')2] ^-p2^p 

^"8C{3)T3J0 wv +(cos26»I'='-6'' ,+0" ,)21[xP , + (cos 251=' -6'' ,Y] 1+eP/T- U^yj 

L -j-g' V u re' re' re' ^' re' re' 
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These equations are coupled differential equations because aF^^i and oF^^, depend on both L^^ 
and Li,^. Recall that = (2 ± z)/2. From Eq.(]93|) the z parameter, which is related to the 
number of sterile neutrinos produced, is governed by: 



^ = 1(1 
dt 4^ 



r^^sin2 2^o-' ^ r^^sin2 2^o-' 



[xl, + (cos 2ef - If,, + al,y] \xl, + (cos W,^' - + 



(130) 



and the evolution of z is governed by an equation similar to the above (but with dP —dP\ 
The above equations can be integrated numerically (following the proceedure mentioned in 
section V). Doing this, we can find the region of parameter space where the L^^"^ asymmetry 
created by the — z/g oscillations does not get destroyed by the — v', oscillations. We 
will solve Eq. ( |129| ) and Eq. ( |130| ) under the assumption that sin^ 26*0^ ~ 1 (i.e. the v^. — v', 
oscillations are approximately maximal). Performing the necessary numerical work, we find 
that L^^^ is created by Vr — v'^. oscillations and not subsequently destroyed by Ve — '^^e oscillations 
for the region of parameter space shown in Figure 5. For definiteness we have taken two 
illustrative choices for sin 2^5^', sin2 205^' = 10-^ 10"^ Note that in our numerical work, we 
have studied the region 10~^ ~ \bvp?^,\l eV''^ ~ 10^. Of course, there will be parameter space 
outside this region where the V^^^ created by v^. — v', oscillations is not destroyed by — v', 
oscillations. However, one should keep in mind that there is a rather stringent cosmology 
bound, vfiy^ ~ 40 e\^[Q (which implies that |5m^g, | ~ 1600 eV'^\ This bound assumes that 
the neutrino is approximately stable, which is expected given the standard model interactions. 
Of course, if there are new interactions beyond the standard model, then it is possible to evade 



this cosmology bound [4£ 



Observe that the region of parameter space where L^^^ survives is somewhat larger than 
our analytical estimate Eq. ( P-28|) . This is partly because the point where L^^ is created occurs 
at a significantly higher temperature than the analytical estimate (see section V for some 
discussion about this point). Note that the quantity \/ /S.ny/ny a T^H"^ /T^,^ ~ T^is- Thus, 
the result that the lepton number is created at a higher temperature than our analytic estimate 
can easily lead to a significant increase in the parameter space. Also, for large sin 261" , the 
magnitude of Ly^ created by — z/g oscillations is considerably larger before the growth of Ly^ 
is cut off by the non-linearity of the differential equation governing its evolution (compare the 
solid line with the dashed or dashed-dotted lines in Figures 3 or Figure 4). Recall that our 
analytical estimate assumed that the creation of L^^^ due to Ur ~ ^'e oscillations had already 
passed the rapid exponential growth phase at the point where the destruction of L'^^'^ due to 
Ve ~ ^'e oscillations reached a maximum. While this latter assumption is generally true for 
small values of sin^ 26*^*^ , it is not true for larger values. In this case, the rate of change of T^l^ 
due to v^. — z/g oscillations will be much larger than our analytical estimate. Consequently, 
the allowed region of parameter space is increased. Thus the result that the allowed region 
of parameter space for sin 2^™ = lO-'^ is significantly larger than the allowed region for 
sin^ 26'J'^ = 10~^ is not unexpected. 

Having established the condition that the — z/g oscillations do not destroy the L'^^^ which 
is created by the — v', oscillations (or — v', oscillations), we must also check that the 
magnitude of L^'^'^ is large enough to invalidate the bound in Eq. (|99D . 



34 



For Sml^, in the range l^m^g, | ~ 10~^ eV'^, the bound Eg. (PU]) arises by requiring that the 
i>e — i/g oscillations do not bring the z/^ sterile neutrino into equilibrium above the kinetic de- 
coupling temperature (~ 3 MeV). The sterile neutrino u'^ will not be brought into equilibrium 
provided that the rate of u'^ production is approximately less than the expansion rate H, i.e. 

r(z/e - u:)/H ^ Ir,^ sin^ 26::' /H ~ 1, (131) 

where we have used Eq.(^ with {siv? t /2Lm) — 1/2|Q. Recall that we are primarily inter- 
ested in the region 1 MeV ~ T ~ 100 MeV, where H ~ h.hT^jMp. Using Eq.(|l9D with a ~ 0, 
the above equation can be re-written in the form 

ye F y < (^^32) 



22 [bl, + 1] 

where we have assumed large mixing i.e. cos 29^" < 1. Recall that b^e' can be obtained from 
Eq.(18). Obtaining the maximum of the left-hand side of Eq. (|132|) leads approximately to 
the bound \5m''\ ~ 10"^ eV. 

In the case where Li,^ is created by — u'^ oscillations, the situation is very different. The 
lepton number L^^ is created at the temperature when bre' ~ 1 (assuming that cos 26^^ ~ 1). 
Denoting this temperature by T^'^ , then as per Eq. ( 




MeV. (133) 

The evolution of this system can be divided into two regions, the region before lepton number 
creation (i.e. T > TJ"^), and the region after the lepton number creation (i.e. T < T^'^). 
In the region before the lepton number is created, a^e' — and Eq. (|132|) holds. We will 
obviously be interested in the parameter space where ^rri^^, is sufficiently large (recall that 
(5m^g/ is related to by Eq.( p.33| ) above) so that L^^ is created at some point above the 



kinetic decouphng temperature T^ec — 3 MeV, of z/g. Let us assume that |5m^g, | is large 
enough so that b^^, ^ 1 for temperatures T > T^^ (which corresponds approximately to, 
|5m^g/| > |(5mgg/|). In this case, 6gg/ + 1 ~ 6gg/ and Eq.( |132| ) can be re- written in the form 

^9 > U.16ml,M^ \ ' sin' 2efy,Gj.Mp 

^ V Q-SV2GfA, 22 ' ^ ' 



where we have used Eq.(^8]) with ~ T /4.I. Observe that the most stringent condition 
occurs for T = T^'^ . Thus taking T = T^^ , using Eq.( |133D , we find that 



\6m. 



2 

re' I > 



15 



|5m; 



3V-2 



4 



(135) 



Assuming maximal Ve — v'e oscillations (i.e. sin 291^ ~ 1) and assummg 
Eq. ( |135| ) imphes that 

|(5m^g,| ~ 10"=^ eV'. (136) 
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Thus provided that this constraint is satisfied, the sterile neutrino, z/g 
equihbrium for temperatures greater than the temperature where L^^ is created, TJ"^ . 

We now need to check that the lepton number created is sufficient to suppress — v'^ 
oscillations for temperatures less than TJ"^ . Demanding that the interactions do not bring the 
sterile neutrino into equilibrium with the muon neutrino, that is again imposing the inequality 
Eq.( |131| ), but this time for T < T^^ where there is significant creation of L^^^ [|ll|], we find 
that 

where the — (+) signs correspond to z/g — z^g ( z/g — z/g ) oscillations. Note that Eq. (|137| ) is only 
required to be satisfied for T > T^^gc — 3 MeV (since we only need to require that the sterile 
neutrinos do not come into equilibrium before kinetic decoupling of the electron neutrinos 
occurs). Once L^^^ is created at T = TJ^' ( where hre! = cos 2^5^ ~ 1), its magnitude will rise 
according to the constraint a^g' ~ 1 (assuming for definiteness that L*^^) > 0). Note that the 
quantities bee',a,ee' are related to foT-g/, follows: 



6ge' A^mL 1 5ml 



re' 



(138) 



After the initial resonance are' ~ 1 while fore' < 1 (and quickly becomes much less than one) 
Thus very soon after the resonance, are' ^ &re' and hence from Eq.( |138| ), a^e' ^ &ee'- As 



before the most stringent bound occurs when T ^ T^'^ , and Eq. ( |137|) leads to approximately 
the same bound as before [i.e. Eq. (|136|) ], since at the point T = TJ"*^ , are' ~ Ke'- 

Finally, we need to check that the oscillations of the Ue, v'^ neutrinos do not significantly 
deplete the number of electron neutrinos for the temperature range, 

0.7 MeV ~ T ~ Tdgg ~ 3 MeV. (139) 

Neutrino oscillations in this temperature range can affect BBN because they will deplete 
electron neutrinos (and anti-neutrinos) and thus modify the temperature when the neu- 
tron/proton ratio freezes out. This effect is generally small unless sin^2^m ~ 10~^P, §]. 
If we demand that sin^2^m ~ 10~^ for this temperature range, then from Eq.(|l9D we require 
\a\ ~ 10 (for the most stringent case of maximal mixing) for this temperature range (or 
\5m?\ ~ 10^^ eV"^). Thus, from Eq.(|18D, \a\ ~ 10 implies 

Recall that for temperatures T ^ Tj (where Tj is the temperature where the change in a due 
to the expansion is larger in magnitude to the change in a due to oscillations, see the earlier 
comments around Eg . (^9]) for some discussion about this), L^^ is created such that are' — 
from this it follows that 



LW=.V=.2xlO- 1^!?^ ^ . (141) 
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Combining Eg. ( 140 ) and Eq.( |141| ), sufficient lepton number will be generated to suppress the 
oscillations in the temperature range Eq. (|139|) provided that 



\SmU ^ lO~'\SmU (^] . (142) 



f 



Note that the temperature Tj is generally less than about 4 MeV [see Eq. ([39| ) for a discussion 
about this]. Thus, Eq.( |142| ) will be easily satisfied given the condition Eq. (|128|) . 



In summary, a consequence of the creation of Lj,^ by u^- — v'^ oscillations is that the large 
angle or maximal — ^'e oscillations will not significantly modify BBN provided that L^^^ does 
not get destroyed by i'e~i-''e oscillations (see Figure 5 for some of this region of parameter space) 
and the condition Eq. (|136|) holds. Thus, it is clear that the oscillation generated neutrino 
asymmetry can weaken the rather stringent BBN bound {\6ml^,\ ~ 10~^ eV'^ for maximal 
mixing) by many orders of magnitude. A consequence of this is that the maximal ordinary- 
sterile neutrino oscillation solution to the solar neutrino problem does not significantly modify 
BBN for a large range of parameters. 

While we have focussed on a particular scenario, our analysis will be relevant to other 
models with sterile neutrinos. For example, assume that there is a sterile neutrino which 
mixes with parameters corresponding to the large angle MSW solution to the solar neutrino 
problem, that is ~ 10~^ eV'^ and sin^ 26'o ~ 0.7[^. This scenario has been "ruled out" 



(assuming negligible lepton number asymmetry) in Refs. ^ . However, if the sterile neutrino 
also mixes slightly with the mu and/or tau neutrino (and such mixing would be expected), 
then these BBN bounds can be evaded provided that |5m^g, | and/or |5m^e'l ~ 0.1 — 1 eV'^. 
Note that the evidence for z/^ — Ue oscillations found by the LSND collaboration suggests that 
|5m^g| ~ 0.3 e^^P, If this is the case then the large angle MSW solution will not lead 
to a significant modification to BBN for a large range of values for sin^ 26q^ . 

We now discuss the possibility that the atmospheric neutrino anomaly is due to large 
angle or maximal muon neutrino - sterile neutrino oscillations. Here, we will denote the 
sterile neutrino by u'^ (this neutrino is expected to be distinct from u'J. Note that the 
possibility that the atmospheric neutrino anomaly is due to large angle or maximal z/^ — z/^ 
oscillations can be well motivated. For example, the exact parity model [|1^ predicts that 
all three ordinary neutrinos mix maximally with mirror neutrinos if neutrinos have mass. 
[See also Ref. [^ for some other interesting models which can solve the atmospheric neutrino 



anomaly through maximal ordinary - sterile neutrino oscillations] . The deficit of atmospheric 
muon neutrinos can be explained if there are — u'^ oscillations with sin 

2^0 ~ 0.5 and 

10-3 ^ \5ml^,\/eV'^ ~ lO^^g, 0. The best fit occurs for 5ml^, ~ 10"^ eV^ and sin^ 2^o - 1 
||2|. However, this parameter range is naively inconsistent with BBN [see Eq.(|l|)] if the lepton 
number asymmetries are neglected. Can the generation of lepton number by ordinary-sterile 
neutrino oscillations reconcile this solution to the atmospheric neutrino anomaly with BBN? 

To study this issue, consider the system consisting of z/,-, z/^, z/^. This system is similar 
to the z/t-, z/e, z/g system that we have discussed above. Doing a similar analysis to the above 
(i.e. replacing v^. and u'^ by z/^ and z/^), we find that the L^^"* asymmetry created by Ur — ^'^ 
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oscillations will not be destroyed by — z/^ oscillations provided that 



(143) 



where A and An^/ni, are given by equations similar to Eq.( |114| ) and Eq.( p.26| ) except that the 
replacements Ue — > y^j., Ae — >• have to be made. Thus, evaluating the resulting expressions 
for A and An^/n^, we find 

5x10- f^f. (144) 



As before, we have made a more accurate numerical study of this problem. If we solve the 
system of equations Eq. ( p-29|) and Eq. ( |130| ), with the replacements z/g, z/g ^/j,,^'^, then we 
can obtain the region of parameter space where the L^^^ created by z/^ — u'^^ oscillations does 
not get destroyed by z/^ — u'^ oscillations. We show some of this parameter space in Figure 6. 
If we assume the best fit of the atmospheric neutrino data, then |5m^^/| ~ 10~^ eV"^ and 

,2 o/gMK 



sin 26q^ ~ 1. Numerically solving the Eqs. ( |129| ) and Eq.( |130| ) (with the replacement of z/g, z^e 
with z/^, z/|j) assuming the best fit parameters, ^ 10^^ eV"^ and sin^ 26'o^ ^ 1, we 

again obtain the region of parameter space where the L^^^ asymmetry is created by z/^ — z/^ 
oscillations and does not get destroyed subsequently by z/^ — u'^ oscillations. Our results are 
shown in Figure 7. As the figure shows, the asymmetry L^^^ created by Ur — u'^ oscillations 
will not be destroyed by z/^ — v'^ oscillations provided that Sm^^, is quite large, i.e., 

\Sml^,\ ^ 30 eV^. (145) 

Recall that our analysis neglects the possible effects of z/^ — z/^ oscillations. It may be possible 
that smaller Sm'^^, are allowed if the Uj. — mixing parameters are large enough. 

The requirement that — v'^ oscillations do not produce too many sterile states implies 

an upper limit on sin^ 2^o^ [see Eq. (|96|) j. This upper limit has been shown in the Figure 
(dashed-dotted line). Also shown in Figure 7 (dashed line) is the cosmological energy density 
bound \5ml^,\ ~ 1600 e1/2[||. 



Recall that the differential equations, Eq. (|129|) are only valid provided that Eq.(|75D holds. 
[We also require Eq.(^) to hold for 5m? = Sm^^,, which is clearly valid for the region of 
parameter space studied]. Note that in our numerical work we found that the condition 
Eq. ([TSp was approximately valid for the points in the allowed region of the figures except for 
the region with relatively large values of sin^ 26'o ~ 10^^. It would be a useful exercise to 
check our analysis by preforming a more accurate study using the density matrix equations, 
Eq.(^6D, modified to incorporate the neutrino momentum distribution. 

The sin^ 29^'^ dependence shown in Figures 6,7 can be understood qualitatively as follows. 
For small sin^ 29^'^ (~ 10~®), the creation of L^'^^ is sluggish which has the effect of delaying 
the point where the destruction of L^'^^ by t'^ — z/^ oscillations reaches its maximum rate. As 
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mentioned earlier, for lower temperatures the rate at which — u'^ oscillations destroy L^'^) 
increases, which has the effect of reducing the allowed parameter space. For larger values of 
sm 

2C (~ 10-8), the maximum rate at which the — v'^ oscillations destroy L*^'^^ occurs 
during the time when L^^"* is still growing exponentially. In this case the system moves rapidly 
away from the — v'^ resonance region. Consequently, the allowed region of parameter space 
is significantly increased. 

Observe that from Eq.( p.35| ), this lepton number will easily be sufficiently large and created 
early enough to prevent the v'^ sterile neutrino from coming into equilibrium given Eq. (|145|) . 
Thus, the large angle or maximal muon - sterile neutrino oscillation solution to the atmospheric 
neutrino anomaly is in fact consistent with BBN for a significant range of parameters. Note 
that the condition Eq. ( |145D implies quite large tau neutrino masses, m^^ ~ 6 eV . Note that if 
the neutrinos are approximately stable (which would be expected unless some new interactions 
exist [^) then there is a stringent cosmology bound of rriy^ ~ 40 el^[0. Although this 
parameter space is not so big, it can be well motivated from the point of view of dark matter 
(since stable tau neutrinos with masses in the range 6 eV ~ mj,^ ~ 40 eV could provide a 
significant fraction of the matter in the Universe). 

If we add the v'^ sterile neutrino to the Vr-, J^e, T^'e system we considered earlier (in connection 
to the large angle ordinary - sterile neutrino oscillation solution to the solar neutrino problem), 
then Ur — u'^ oscillations will also generate L^^^ in a similar manner to the way in which z/^ — z/g 
oscillations generated L'^^\ Consequently, the bounds on (5m^g,, sin^26'Q'^', can alternatively 
be considered as bounds on Sm"^^, , sin^ 29q^ . Of course, we only need to require that either 
Sml^,, siii^ 29q^' or sin^26'o'^ satisfy the bounds derived. Similarly, we can add the v'^ 

sterile neutrino to the lyr^^^i, ^'^ system and analogous reasoning leads to the conclusion that 
the bounds on dm^^,, sin^ 26*0^^' can alternatively be considered as bounds on Srn^^/, sin^ 26q^ . 

Observe that with 6mz^,, sin 2^5^' or Sml^,, sm^2ei'' in the range identified in Figure 7 
(where the atmospheric neutrino anomaly is explained by large angle i/^ — z/^ oscillations 
without significantly modifying BBN) the solar neutrino problem can also be solved for the 
entire parameter space [Eq.(p8D], without significantly modifying BBN. Alternatively one can 
argue that the present data may allow the u'^ to come into equilibrium with the ordinary 
neutrinos and still be consistent with BBN and thus we only require the less stringent 
bounds given in Figure 5. Clearly this is a possibility at the moment. Note however, that 
for the case of the exact parity symmetric model[T^, where the mirror neutrinos interact 
with themselves, this way out is not possible. This is because if the mirror muon neutrino 
is brought into equilibrium above the kinetic decoupling temperature (which is about 5 MeV 
for muon neutrinos) then the mirror weak interactions will bring all three mirror neutrinos 
together with the mirror photon and mirror electron-positron into equilibrium (which would 
lead to about 9 effective neutrino degrees of freedom during nucleosynthesis). For the case 
of mirror neutrinos it seems to be necessary to ensure that the mirror muon neutrino is not 
brought into equilibrium in the first place. 

Note that in our previous analysis, we have assumed that the sterile neutrino is truly 
sterile and does not interact with the background. In the special case of mirror neutrinos. 
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the mirror neutrinos are expected to interact with the background because they interact with 
themselves [0]. In general the effective potential describing coherent forward scattering of the 
neutrino with the background has the form V = Va — V^. For truly sterile neutrinos, = 
(as has been assumed hitherto). For mirror neutrinos is non-zero. Denoting the mirror 
neutrinos by z/^, then for the case of z/^ — z/^ oscillations we will denote the effective potential 
by 

V = V^- (146) 

where Va is given by Eq. (|l^ and is the effective potential due to the interactions of the 
mirror neutrinos with the background. The mirror effective potential can be expressed in 
an analogous way to Va, that is there is a part which is proportional to mirror lepton number 
and a part which is independent of mirror lepton number, 

V; = i-a'^ + b'nK (147) 

If the number of mirror neutrinos is much less than the number of ordinary neutrinos then 
6' ~ 0. [Note that the 6-part of the effective potential is proportional to the number densities 
of the background particles. This dependence is not given in Eq.([T5|) since for this equation 
the number densities were set equal to their equilibrium values]. The parameter a' has the 
form 

a".Z^%^, (148) 

where L'^^^ is given 

L'(^) = L,. + L,, + L,. + L,, + r]\ (149) 

where L^i^ are are the mirror lepton numbers, which are defined by L^/^ = {n^/^ — np'^)/n^ 
(note that is the number density of ordinary photons) and rj' is a function of the mirror 
baryon/electron number asymmetries [which is defined analogous to Eq.(0)]. We will assume 
that 1]' is small and can be approximately neglected. Since ordinary + mirror lepton number 
is conserved (and we will assume that it is zero), it follows that 

L^^ + L^^ + L^^ + L^, + L^. +L^,=0. (150) 

From the above equation, it follows that a' is expected to be of the same order of magnitude 
as a. In the case of the z/,-, z/^, u'^ system, the effect of the mirror - neutrino effective potential 
can be accounted for by simply replacing L*^^'"^) in Va by 

Lir) ^ L^r) _ ^/(M) ^ 2L,^ + L,^ - ~ AL,^ + 3L,^, (151) 

where we have used Lj^^ ~ Lj^/^ ~ and Eq.( |150| ). Thus, from the above equation, assuming 
that negligible L^^^ is produced, we see that \a^^i\ ~ |i?|a,-^/| (where R = |5m^^,/5m^^, |). The 
factor of 3/4 replaces the factor of 1/2 that we had earlier (for the case where u'^ or z/^ were 
sterile neutrinos). This difference will increase the region of allowed parameter space, because 
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it will make T/^^' closer to the point where L^^ is initially created. At this point dLy^/dt can 
be significantly enhanced because it is very close to the resonance (also note that dL„^/dt 
will be suppressed because it is proportional to 1/T^). 

Finally, observe that another important feature of mirror neutrinos is that the mirror inter- 
actions can potentially bring all three of the mirror neutrinos into equilibrium with themselves 
as well as the mirror photon and mirror electron positron. [However the temperature of the 
mirror particles will generally be less than the temperature of the ordinary particles if the os- 
cillations satisfy Eq.(p6|)]. Detailed studies involving mirror neutrinos will need to incorporate 
this. We leave a more detailed study of mirror neutrinos to the future . 



VII. Conclusion 

In summary, we have studied the phenomenon of neutrino oscillation generated lepton 
number asymmetries in the early Universe in detail. This extended study clarifies the origin 
of the approximations adopted in the earlier work||12||. We have also studied the effects of 



the thermal distribution of the neutrino momenta and non-negligible sterile neutrino number 
densities. 

In the unrealistic case where the neutrino momentum distribution is neglected, the evolu- 
tion of L,y^ can be approximately described by seven coupled differential equations [Eqs. (p6D ], 
which can be obtained from the density matrix. We showed in section III that these equations 
can be reduced to a single integro-differential equation (we show in the appendix that the same 
equation can be obtained from the Hamiltonian formalism). In general, the density matrix 
equations cannot be solved analytically, and must be solved numerically. However, if the 
system is sufficiently smooth (the static limit), then the integro-differential equation reduces 
to a relatively simple first order ordinary differential equation [Eq.(p3D]. This equation gives 
quite a reasonable description of the evolution of L^,^ , except possibly at the initial resonance 
where significant generation of L^,^ occurs. We show that when the thermal distribution of 
the neutrino momenta is incorporated several important effects occur. One of these effects 
is that the creation of lepton number is much smoother. This allows a considerable compu- 
tational simplification, because it means that the static approximation can be a reasonably 
good approximation, even at the resonance for a much larger range of parameters. This means 
that L,^^ can be accurately described by the relatively simple first order differential equation 
(modified to incorporate the neutrino momentum distribution). This equation is given by 
Eq. (p4D , expressed as a function of the number distribution of sterile states. In section V, we 
showed that the number distribution of sterile neutrino states satisfied a first order differential 
equation [ Eq.(^)] which must be integrated for each momentum step. 

We first applied our analysis to obtain the region of parameter space where large neutrino 
asymmetries are generated. This region of parameter space is given in Eq. (p7|) . This analysis 
included the effects of the neutrino momentum distribution which was neglected in earlier 
studies P^, |TB[. We also examined the implications of lepton number generation for the BBN 



bounds for (5m^,sin^2^o for ordinary-sterile neutrino mixing. There are two ways in which 
the creation of lepton number can modify the BBN bounds. One way is where the Ua — 
Ps oscillations themselves produce Ly^ thereby suppressing the number of sterile neutrinos 
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produced from the same oscillations. The other way is where the — fs oscillations create 
L^i^ which thereby suppresses production from z/q, — oscillations. The bound for the 
former case is given in Eg. (PBD , while the latter case studied in section VI, in the context of 
the maximal vacuum oscillation solutions to the solar and atmospheric neutrino problems. The 
maximal vacuum oscillation solution of the solar neutrino problem assumes that the electron 
neutrino is approximately maximally mixed with a sterile neutrino. For a large range of 
parameter space, the maximal mixing leads to an energy independent factor of two reduction 
in the solar neutrino fluxes. This leads to a reasonably simple predictive solution to the solar 
neutrino problem which is supported by the experiments. However, most of the parameter 
space for this solution is inconsistent with standard big bang nucleosynthesis (BBN) if the 
lepton numbers are assumed to be negligible |^, ^ ||. We showed that there is a large region 
of parameter space where the oscillations generated lepton number in such a way so as to allow 
the maximal vacuum oscillation solution to the solar neutrino problem to be solved without 
significantly modifying BBN. The allowed parameter space is given in Figure 5. We also 
showed that there is a range of parameters where the lepton number is generated so that the 
large angle muon - sterile neutrino oscillation solution to the atmospheric neutrino anomaly 
does not lead to any significant modification of BBN. This parameter space is illustrated in 
Figure 6 and Figure 7. 

We finish with a speculation. One of the mysteries of cosmology is the origin of the 
observed baryon asymmetry of the early Universe. In principle, it may be possible that the 
baryon asymmetry arises from a lepton number asymmetry. The lepton number asymmetry 
can be converted into a baryon number asymmetry through sphaleron transitions at or above 
the weak phase transition. It may be possible that a small lepton number asymmetry arises 
from the mechanism of ordinary-sterile neutrino oscillations, which is seeded by statistical 
fluctuations of the background. One interesting feature of this possibility is that the baryon 
number asymmetry would not be related to the CP asymmetry of the Lagrangian. Instead 
the origin of matter over anti-matter would be due to a statistical fluctuation which is then 
amplified by neutrino oscillations. However before this speculation can be checked, it would 
be necessary to work out the effective potential at high temperatures (T ~ 250 GeV) and 
study the phase transition region. 
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Appendix 

The purpose of this appendix is to show that Eg . (|60|) can be derived from the Hamiltonian 
formalism. In applying this formalism, we will assume that the rate at which collisions 
collapse the wavefunction (i.e. the rate of measurement of whether the state is a weak or 
sterile eigenstate) is given by the damping frequency which is half of the collision frequency. 
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For further discussion of this point see section II and Ref. [ 

The expectation value that an initial weak-eigenstate neutrino z/^ has oscillated into a 
sterile state after r seconds will be denoted by \ip'g{t, r) p (where t is the age of the Universe). 
The average probability that an initial weak eigenstate has oscillated into a sterile state can 
be obtained by averaging the quantity |?/'^(t,r)p over all possible times r (weighted by the 
probability that the neutrino has survived r seconds since its last "measurement"). This 
average has the form 

1 



uq Jo 



(152) 



where cuq is the mean time between measurements. According to Ref.[^], ujq = 1/D = 

If we denote the analogous quantity for anti-neutrinos by (| ip^ (t, r)p), then the rate of 
change of lepton number can be expressed as 



dt 



8^ ^ " 2 



3 d{n{t)) 
8 dt 



(153) 



where 



{m) = mt)\')-{\i^sm- 

Note that the first term in Eq. (|153|) represents the rate of change of lepton number due to 
collisions (which produce sterile neutrino states). The second term represents the rate of 
change of lepton number due to the oscillations between collisions. 

In the adiabatic limit, the transformation ~^ and Lq — > Lm diagonalizes the Hamil- 
tonian. In this limit, the mean probability (|^s(t)P) is given by: 



(|V';(t)P)=sin2 2^^(sin^ 



T 



2Lr 



(155) 



Note that in the static limit, d{Vt{t)) / dt = and hence Eg . (pOD results. However, in the 
expanding Universe which is non-static, the above equation is not generally valid (although 
it turns out that it is a good approximation for oscillations away from resonance where the 
system changes sufficiently slowly and even at some resonance regions which are sufficiently 
smooth). To calculate the probability (|V's(^)P) in the general case, we go back to the funda- 
mental Hamiltonian equations 



dt [ ^; 



1 

2^ 



where 



and 




6m' +4P 



/ cos^o sin^o \ ,Tj\ ^ {b±a)Sm' 
^ [ -sin^o cos^o ; ' ^ ' 2p 



(156) 



(157) 
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where the — (+) sign corresponds to neutrino (anti- neutrino) oscillations. Expanding out 
Eq.(|l5|), we find: 



.#0 



dt 



where 



a = ^(26 ± 2a -cos 2^0), P 



sin26'o, 1 = —. — cos 26^0. 



(159) 



(160) 



4p 2p Ap 

If we divide the equations Eqs.( |1590 by ipa and ip's respectively, then they can be combined 
into the single differential equation: 



where = ip'^/ipa and 



A = 7 — a 



6m^ 



[cos 26q — 6 ± a). 



(161) 



(162) 



The +(— ) sign in the above equation corresponds to — {^a — Vg) oscillations. Note that 
Ws\^ = i+\ki\-2 ■ If non-linear term (\1/^) can be neglected, then the solution for constant 



a, /3, 7 is 



VI/ (t) 



2A 



,-i(A)(t-i*) 



(163) 



with boundary condition \l'(t*) = 0. Introducing the variable r = t — t*, and evaluating 
|^(t,r)|2 we find 

|vl/(t,r)p = ^sin^[^], (164) 

which is approximately sin^ 2^^ sin^ r/2Lm provided that |\E'p ^ 1. 

In the general case where a, (3 and 7 are not constant, the general solution is (where we 
have again neglected the non-linear term): 



where 



A(t') 



xdtr 



(165) 



(166) 



and the boundary condition \E'(t*) = has again been taken. One may easily verify that 
Eq.( |165| ) is indeed the solution by directly substituting it into Eq.( |161| ). The probability that 
a weak-eigenstate aXt = t* has oscillated into a sterile eigenstate at time t is thus 



^(t)|2~ ^ r r COS \ r xdt' 

4 Jt* Jt* VJt2 



2 rt 
t* Jt* 



dtidt2 



(167) 



where we have assumed that (3 is approximately constant over the interaction time scale t — t*, 
so that it can be taken outside the integral. This step is a good approximation provided 
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T ~ 2 MeV[^. Again defining the quantity t = t — t* (recall that r is the time between 
measurements), and averaging |\E'(t, r)p over r, with the appropriate weighting factor, we find 
that 



4:Uo Jo 



-t/ujq 



COS 



t-T Jt-T 



Xdt' 



t2 



dtidt2dT. 



Integrating this equation by parts (with respect to the r integration), we find: 



e~^/'^o cos 



Jt-T 



Xdt' 



t~T 



dtidr, 



:i68) 



(169) 



where we have used the fact that e ~ 0[Q. The analogous quantity for anti-neutrinos, 

(I ^ (^)P), can similarly be defined. Recall that the functions (|\E'(t)p), (| ^ (t)p) are related 
to the rate of change of lepton number through Eq. (|153|) : 



dL,. 



dt 



dt 



(170) 



Evaluating d{\^{t)\'^) /dt we find: 

dt 2 



/o*e-M.(cos [lUXdf 



- 1 dr 



+ f /o e--/-"A(i_.) Jl^ sin Xdt'] dt.dr, (171) 
where we use the notation that X(^t-T) denotes A evaluated at the point {t — t). Dividing 



Eq. ( |169| ) by and integrating Eq. 
we find 



by parts we find (with respect to the r integration), 

— (I^WH = C Te-"/"" fl - X^t-r) f sin [ r Xdt] dt,) dr. (172) 
ojq 2 Jo V Jt-T \_Jt-T J / 

Adding the above two equations and subtracting the analogous term for anti-neutrinos, we 
obtain the following rather compact expression for the rate of change of lepton number: 



dL, 



dt 



-3^ 
16 Jo 



r/u}o 



COS 



Xdt' 



t-T 



— COS 



Xdt' 



t-T 



dr. 



(173) 



where A is defined similarly to A except that a — *• —a. Note that the total contribution to the 
rate of change of lepton number is in fact simpler than either of the two separate contributions 
coming from collisions and oscillations between collisions. Equation ( |173| ) can be re-written 
(using a trigonometric identity) 



dL^^S^ 
dt 8 



e-^/^° sin 



X+dt' 



t-T 



sm 



/* x-dt" 

Jt-T 



dr. 



(174) 



where A^ = (A ± A)/2. Note that this is exactly the same equation that we derived in section 
III [Eq. (|60|) ] from the density matrix equations. 
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Table Caption 



Summary of the predictions for the chlorine and gaUium experiments assuming 1) standard 
electro-weak theory (i.e. no new physics) 2) that the electron neutrino oscillates maximally 
into a sterile state (maximal mixing model) and 3) the experimental measurements. All 
numbers are in units of SNU. 

Figure Captions 

Figure 1. The evolution of the i/^ — z/^ (or Vj. — Vg) oscillation generated lepton number 
asymmetry, Ly^^ (or L^^). We have taken by way of example, the parameter choices 5m? = 
— 1 eV^^, sin^ 26*0 = 10"^ (sin^ 26^0 = 10~^) for the bottom two curves (top two curves). The 
solid lines represent the results of the numerical integration of the density matrix equations 
[Eq.(^BD], while the dashed lines result from the numerical integration of Eq.(PBD. 

Figure 2. The evolution of the z/^ — (or u^- — Ug) oscillation generated lepton number 
asymmetry, L^^ (or L^^). In this example we have taken the parameter choices, = 
—1000 eV'^, sin^ 26*0 = 10~^ (sin^ 26*0 = 10~®) for the bottom two curves (top two curves). As 
in Figure 1, the solid lines represent the results of the numerical integration of the density 
matrix equations [Eq.(|46|)], while the dashed lines result from the numerical integration of 
Eq.(^. 

Figure 3. The evolution of the z/^ — Ug (or u^- — Ug) oscillation generated lepton number 
asymmetry, L^^ (or L^^). In this example we have taken the parameter choices (5m^ = 
-1 e1/2, sin^ 200 = 10"^ (dashed line), sin^ 2^o = 10"^ (dashed-dotted line) and sin^ 2^o = 10"* 
(solid line). These curves result from integrating the coupled differential equations Eq.(p4|, 
p3D , which in contrast to Figures 1,2, incorporate the momentum distribution of the neutrino. 
They also incorporate the effect of the non-zero number density of the sterile neutrinos which 
are produced by the oscillations. 

Figure 4. Same as Figure 3 except that (5m^ = —1000 eV^,sin^26'o = 10^® (dashed line), 
sin2 2^o = 10"^ (dashed-dotted line) and sin^ 2^0 = 10~^ (sohd line). 

Figure 5. Region of parameter space in the — 5m^g,, |5mgg/|, plane (assuming sin^ 26'q'^' ~ 1) 
where the L^^^ created by i^r — ^'e oscillations does not get destroyed by z/g — z/^ oscillations. The 
solid line corresponds to sin^ 26'™ = 10~^, while the dashed line corresponds to sin^ 26'J'^ = 
10~^. Note that similar results hold for z/^ — z/g oscillations by replacing Ur — ^ z/^. 

Figure 6. Region of parameter space in the —5m^^,, \5mj^^,\ plane where the L^^^ created 
by Ut- — p'^ oscillations does not get destroyed by z/^ — p'^ oscillations (assuming sin^ 26q^ ~ 
1). The solid line corresponds to sin^ 2^0*^' = 10^^, while the dashed line corresponds to 
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sin^ 29 = 10 ^. Note that similar results hold for — v'^ oscillations if both i/^ and z/g exist. 

Figure 7. Region of parameter space (sin^ 26'o'^ , —6ml^,) where the L^^^ created by Vr — ^'^ 
oscillations does not get destroyed by — v'^ oscillations. This region which in the figure is 
denoted by the "Allowed region" is all of the parameter space above the solid line. We have 
assumed that sin^ 2^0^^ ~ 1 and | = 10^^ eV"^ (which is the best fit to the atmospheric 

neutrino data). Also shown (the dashed line) is the cosmology bound m^,^ ~ 40 eV (which 
implies |(5m^^, | ~ 1600 eV"^), which is required if the neutrino is sufficiently long lived. The 
dashed-dotted line is the BBN bound, Eq.(^). 
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Table 1 



Prediction/Expt 


Chlorine 


Gallium 


Standard Electro-weak theory 
Maximal mixing model 
Experiment 


4.5 ±0.5 
3.7 ±0.4 
2.78 ±0.35 


123+« 
6511 
71 ±7 



53 




Temperature (MeV) 




Temperature (MeV) 



O 

-I— > 



10" 



10"^ ---- .figure 3 



10"' 
10"' 
10"^ 
10"^ 



10"^. 



10-^^ 



10 



11 



10 



-12 




10 15 
Temperature (MeV) 



25 




Temperature (MeV) 



